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CHAPTER  I   INTRODUCTION 

The  overall  objective  of  this  study  is  to  carry  out  a  comprehensive 
investigation  of  the  performance  of  the  semif luidlzed-bed  filter.  This 
filter  has  been  proposed  as  a  remedy  for  the  inefficient  operation  of 
deep  bed  filters  (Wen  and  Fan,  1979;  Hsu  and  Fan,  1981;  1982;  Fan  and  Hsu, 
1984).   Because  the  filtration  in  a  semifluidized  bed  involves  the  deep 
bed  filtration,  a  reasonable  model  of  semifluidized-bed  filtration  requires 
the  incorporation  of  a  deep  bed  filtration  model.  Therefore,  a  companion 
study  on  the  deep  bed  filter  has  also  been  carried  out. 

Deep  bed  filters,  also  known  as  granular  filters,  are  extensively 
used  in  purification  of  potable  water  and  in  treatment  of  waste  water 
(see,  e.g.,  Fuller,  1898;  Dunbar  and  Calvert,  1908;  Hall,  1957).   Much 
effort  has  been  spent  to  study  various  types  of  deep  beds  under  a  wide 
range  of  operating  conditions.   It  appears  that  so  far  little  has  been  done 
to  analyze  and  model  mechanistically  the  results  in  a  unified  manner. 

The  use  of  stochastic  models  to  simulate  the  deep  bed  filtration  has 
been  proposed  by  Litwiniszyn  (1963;  1966;  1967;  1968a;  1968b;  1969).  These 
stochastic  models  can  be  viewed  as  the  intermediate  between  the  phenomenolo- 
gical  models  (Iwasaki,  1937;  Ives,  1960;  1961;  Camp,  1964),  in  which  the 
form  of  model  equations  is  assumed  from  prior  knowledge  of  the  phenomena, 
and  the  trajectory  analysis  models  (O'Meliaand  Stum,  1967;  Yao  et  al. , 
1971;  Payatakes,  1973;  Rajagopalan  and  Tien,  1976;  1979;  Pendse  et  al. , 
1978),  In  which  the  trajectories  of  the  particles  are  determined  pains- 
takingly  from  the  force  balance  equations.  The  filtration  process  is 
complicated  and  chaotic  in  nature,  and  the  stochastic  models,  derived 
through  probability  considerations,  often  generate  parameters  that  appear 
to  describe  adequately  the  behavior  of  the  filter. 
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Litwiniszyn  (1963)  has  used  a  pure  birth  process  where  the  number 
of  blocked  pores  in  a  unit  volume  of  the  bed  is  considered  as  the  model 
variable.  Hsu  (1981),  and  Hsu  and  Fan  (1984)  have  employed  the  pure  birth 
process  to  model  the  pressure  drop  through  the  entire  bed;  in  their  model 
the  number  of  blocked  pores  is  related  to  the  pressure  drop  by  assuming 
that  the  scouring  of  deposited  particles  is  negligible. 

One  group  of  research  workers (see,  e.g.,  Mintz ,  1951;  Mintz  et  al. , 
1967)  have  considered  that  deposits  accumulated  in  the  filter  medium  have 
a  structure  with  uneven  strength.   This  structure  is  partially  destroyed 
by  the  hydrodynamic  forces  of  the  flow  of  water  through  the  medium,  which 
increase  with  an  increase  in  the  pressure  drop;  a  portion  of  the  particles 
less  strongly  linked  to  other  particles  is  detached  from  the  grains.  Based 
on  this  observation,  Litwiniszyn  (1956)  has  developed  a  birth-death  model. 
In  the  present  work,  this  birth-death  model,  taking  into  account  both 
blockage  of  pores  by  suspended  particles  and  scouring  of  deposited  particles, 
is  combined  with  Carman-Kozeny  equation  for  the  pressure  drop  through  the 
packed  bed  to  describe  its  temporal  change  during  a  filtration  run. 

The  applicability  of  the  birth-death  model  to  the  experimental  data 
has  amply  been  demonstrated.   It  has  been  found,  however,  that  in  some 
situations  the  pressure  drop  characteristics  of  the  deep  bed  filtration 
deviate  from  the  behavior  predicted  by  this  model.   To  elucidate  this 
deviation,  a  different  stochastic  model,  namely,  the  second-order,  pure 
birth  process,  has  also  been  proposed. 

The  operation  of  deep  bed  filters  is  often  inefficient;  it  is  typically 
characterized  by  the  shallow  penetration  of  suspended  solids  and  rapid 
clogging  in  the  bed.  Fortunately,  these  disadvantages  can  be  overcome  by 
the  use  of  semif luidized-bed  filters.  An  experimental  study  on  the  perfor- 
mance of  the  semif luidized-bed  filter  has  been  carried  out  to  observe  the 
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phenomena  involved  in  the  semifluidized-bed  filtration  and  to  determine 
advantages  achievable.  To  this  end,  this  work  is  divided  into  the  following 
chapters . 

In  Chapter  II,  a  review  of  the  semifluidized  bed  literature  is  presented. 
A  description  is  given  of  the  methods  of  operating  a  semifluidized  bed. 
Numerous  articles  pertaining  to  the  fundamental  aspects  of  semif luidization 
are  reviewed;  mass  balance  equations  and  empirical  correlations  relating 
the  structure  of  the  bed  to  operating  parameters  are  presented.  Also 
described  in  this  chapter  are  some  applications  of  the  semifluidized  bed. 

The  Master  equation  and  one-step  processes  are  described  in  Chapter  III. 
The  Master  equation,  governing  the  probability  balance  in  a  stochastic  process, 
is  relatively  easy  to  manipulate;  it  can  be  related  directly  to  the  physical 
process  to  be  modeled.  The  one-step  processes  are  a  special  class  of 
stochastic  processes;  it  is  characterized  by  the  occurrence  of  birth  and 
death  events.   The  one-step  operator  technique  is  presented  to  obtain  useful 
information,  such  as  the  mean  and  variance  of  the  random  variable,  from  the 
Master  equation  for  a  one-step  process. 

Chapters  IV  and  V  deal  with  the  stochastic  modeling  of  the  deep  bed 
filtration  under  the  constant  flow  condition.   Specifically,  the  pressure 
drop  dynamics  in  the  deep  bed  filter  are  simulated  in  Chapter  IV  by  the  linear, 
birth-death  process.  The  model  is  relatively  simple  in  that  the  entire  bed 
is  assumed  to  be  spatially  lumped.  The  available  experimental  data  from 
numerous  researchers  are  analyzed  and  simulated  in  the  light  of  this  model. 

In  some  situations,  the  pressure  drop  dynamics  in  the  filter  bed  deviate 
from  the  behavior  predicted  by  the  linear,  birth-death  process.   To  explain 
this  deviation,  a  second-order,  pure  birth  model  has  been  proposed  in  Chapter  V. 
The  bed  is  divided  into  several  spatially- lumped  compartments;  each  depicted 
by  one  of  the  present  models. 
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A  large  number  of  experimental  runs  were  conducted  with  deep  bed 
filters;  analysis  and  simulation  of  the  results  based  on  the  present  models, 
namely,  the  linear,  birth-death  and  second-order,  pure  birth  processes,  are 
presented  in  Chapter  VI.  The  primary  concern  in  this  work  is  to  investigate 
how  the  size  distribution  of  suspended  solids  and  suspension  concentration 
affect  the  performance  of  the  deep  bed  filtration. 

In  Chapter  VII,  an  experimental  study  of  the  semif luidized-bed  filter 
is  presented.  The  results  show  that  the  filtration  capacity  of  the  semi- 
fluidized  bed  filter  is  considerably  greater  than  that  of  an  equivalent  deep 
bed  filter.   In  addition  to  pressure  drop  and  effluent  quality  data,  a 
series  of  photographs  is  included  to  describe  and  explain  the  mechanisms 
of  the  semifluidized-bed  filter,  giving  rise  to  its  superior  performance 
over  the  deep  bed  filter. 

In  the  final  chapter,  significant  conclusions  of  this  work  are 
summarized  and  possible  area  of  future  research  is  suggested. 


CHAPTER  II   PHENOMENA  OF  SEMFLUIDIZATION 

The  idea  of  semif luidization  was  conceived  by  L.  T.  Fan  (1958).   Semi- 
fluidlzation  refers  to  the  phenomenon  created  when  a  mass  of  fluidized  solid 
particles  is  partially  restricted  from  above  with  a  porous  restraining  plate, 
giving  rise  to  the  formation  of  a  fixed  bed  and  a  fluidized  bed  in  series 
within  a  single  containing  vessel. 

A  typical  semif luidized  bed,  as  shown  in  Figure  II-l,  consists  of  a 
vertically-extending  column  providing  a  bed  section  for  the  fluidized  or 
fixed  solid  particles,  a  fluid  distributor  plate  for  introducing  the  fluid 
(gas  or  liquid)  from  the  bottom  of  the  bed,  and  a  movable  upper  restraining 
plate  which  passes  the  fluid  but  prevents  the  passage  of  the  solid  particles. 
The  degree  of  fluidization  occurring  in  the  bed  can  range  from  complete 
fluidization  (conventional  fluidized  bed)  to  no  fluidization  (conventional 
fixed  or  packed  bed)  by  lifting  or  lowering  the  upper  restraining  plate, 
while  the  fluid  is  being  introduced  at  a  flow  rate  greater  than  the  minimum 
fluidization  velocity  for  the  bed.  At  a  given  fluid  velocity  greater  than 
the  minimum  fluidization  velocity  and  less  than  the  maximum  fluidization 
velocity,  the  upper  restraining  plate  can_ be  raised  above  a  point  where  the 
entire  bed  is  completely  fluidized.  When  the  restraining  plate  is  lowered, 
the  uppermost  portion  of  the  bed  is  first  fixed,  and  as  the  lowering  of  the 
restraining  plate  continues,  a  progressively  increasing  portion  of  the  bed 
becomes  fixed.   Simultaneously,  the  fluidized  portion  of  the  bed  right  above 
the  distributor  plate  becomes  progressively  smaller  until,  finally,  a 
condition  is  reached  where  the  entire  bed  is  fixed  (Fan  e_t  al.  ,  1959;  Fan 
and  Wen,  1960;  1968).  Thus,  for  a  given  system,  any  degree  of  fluidization 
between  the  limits  of  a  completely  fluidized  state  and  a  completely  fixed 
state  may  be  attained  by  controlling  the  fluid  flow  rate  and/or  the  restrain- 


Ing  plate  position.   The  mobility  of  the  upper  restraining  plate,  however, 
is  not  essential  to  the  formation  of  a  semifluidized  bed. 

In  case  where  the  density  of  the  bed  particles  is  less  than  that  of 
the  fluid,  the  structure  of  the  semifluidized  bed  is  reversed.  The  fluid 
is  now  introduced  from  the  top  of  the  bed.   The  force  of  the  fluid  on  the 
particles  counteracts  the  buoyancy  of  the  particles  in  maintaining  fluldi- 
zation.   The  fixed  section  of  the  bed  lies  below  the  fluidized  section  and 
adjacently  above  the  restraining  plate  as  shown  in  Figure  II-2. 

STRUCTURE 

A  simple  material  balance  yields  a  relationship  expressing  the 
height  of  the  fixed  or  packed  section  in  a  semifluidized  bed  as  a  function 
of  operating  conditions  (Fan  et  al. ,  1959;  1960;  1961).  Under  the  assumpr- 
tions  that  the  particles  are  uniformly  distributed  in  a  fluidized  bed,  the 
movement  of  a  particle  in  the  suspension  is  completely  independent  of  other 
particles,  and  the  suppression  of  the  expansion  of  the  suspension  and 
subsequent  formation  of  the  packed  section  does  not  change  the  average 
particle  distance  in  the  fluidized  section;  and  the  voidage  of  the  packed 
section  is  constant  and  is  that  of  the  least  dense  static  bed  under  resting 
conditions,  the  following  equation  for  the  height  of  the  packed  section  can 
be  obtained  as  (referring  to  Figure  II-l) 


h  p  (1  -  c   )  -  [h 
pa  s       pa       pa 


^a  -  (''f  -  ^>  77^^  "I-l> 


where  h   is  the  height  of  the  packed  section,  p  ,  the  density  of  the  solid 

particles,  £   ,  the  porosity  of  the  packed  section,  h  ,  the  height  of  the  fully 
pa  t 


fluidized  bed   ,  h,  the  height  of  the  semif luidized  bed,  and  g^,  the 
porosity  of  the  fluidized  section.   It  follows  then,  the  weight  fraction 
of  the  solid  particles  in  the  packed  section,  X,  can  be  expressed  as 

P  A(l  -  c   )(h,  -  h)(l  -  €,) 
y  .   s       pa    f f_  (II-2) 
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where  A  is  the  cross-sectional  area  of  the  bed  and  W,  the  weight  of  the  bed 
particles.  The  observed  heights  of  the  packed  section  have  been  plotted 
for  comparison  with  the  values  calculated  in  accordance  with  Eq.  II-l  by 
Fan  and  Wen  (1961)  as  shown  in  Figure  II-3. 

FLOW  EQUATION 

Dimensional  analysis  of  the  continuity  and  moranentum  equations  for 
particulate  fluidization  led  to  the  expression  (Beranek,  1958) 


—  )  "  0  ,  .        (II-3) 


where  e  is  the  porosity  of  the  static  bed  under  the  least  dense  condition, 

and  u,  u  ,  and  u  are  the  superficial  fluid  velocity,  the  minimum  fluidization 
mf      t 

velocity  in  terms  of  the  superficial  velocity  and  terminal  velocity  of  the 
particles,  respectively.  Fan  and  Wen  (1961)  have  applied  similar  analysis 
to  semifluidization  and  obtained  the  relationship 

h  -  h     C  -  C 

pa    c    or 

where  G       and  G     are  the  minimum  fluidization  velocity  and    the   free    fall 
mf  t 

terminal  velocity,  respectively,  in  mass  flow  rate  per  unit  area  of  the 
empty  bed.  Experimental  data  from  solid-liquid  (Fan  and  Wen,  1961)  and 
solid-gas    (Wen  £t  al. ,    1963)    systems   substantiate   this    relationship.     A 
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two  parameters  fit  by  Kurian  and  Rao  (1970)  have  yielded  an  explicit 
form  of  Eq.  II-4  as 

G-G,          h-h,.  '■      ,, 

[-— ^]  -0.61  (-—HI) -1-2  (11-5) 


On  the  other  hand,  Rac  and  Sarma  (1973)  introduced  the  minimum  semifluidi- 

zation,  G    ,  in  place  of  the  minimum  fluidization  velocity  in  Eq .  II-4, 
osf 

and  developed  the  following  one  parameter  expression; 

'"  "  ^     p*^  '  '^osf.0.2  ■  (11-6) 


Equations  II-4  through  II-6  define  the  state  of  semif luidization. 

Other  investigators  (Babu  Rao  and  Doraiswamy,  1967;  Roy  and  Sarma,  1973; 

Roy  and  Sen  Gupta,  1975;  Roy,  1971)  have  correlated  empirically  the  state 

of  semif luidization  in  various  forms  of  the  function 

h     D    0       h 

c    o    p    f       c 

The  correlations  are  listed  in  Table  ll-l. 

The  minimum  semif luidization  velocity,  defined  as  the  superficial 

fluid  velocity  at  which  the  formation  of  the  packed  section  initiates, 

can  be  obtained  from  Eq.  II-5  or  11-7  by  allowing  h  -»0  .   Poddar  and  Dutt 

pa 

(1969)  have  employed  the  voidage  function  of  Wen  and  Yu  (1966)  and  the 
relationship  between  e   and  the  voidage  function  (Wilhelm  and  Kwauk  1948; 

Lewis  et  al. ,  1949;  Richardson  and  Zaki,  1954)  to  obtain  the  condition  of 

minimum  semif luidization  given  below. 
1.687 
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the  Galileo  number,  and  (jj  ,  the  sphericity  of  particles.  Empirical 
correlations  are  also  available  in  the  form  of  (Roy  and  Sarma,  1972; 
Roy,  1971;  Roy  and  Sen  Gupta,  1975) 

G   .  G   ,  D  p 

f[^  or  -^  ,  /  ,  A  R]  -  0                     (11-9) 

G  ,  G   .  d  p^ 

ml  msf  p  f 


Table  II-2.  The  minimum  fluidization  velocity,  G  ,,  in  Roy  and  Sarma 's 

mt 

correlations  can  be  calculated  from  Leva's  (1959)  simplified  equation  in 
FPS  units: 

=.f  '  "«  ''p  0.88 "^-^°) 


The  maximum  semif luidization  velocity,  defined  as  the  superficial 
fluid  velocity  at  which  all  solid  particles  are  supported  by  the  fluid  in 
the  packed  portion  of  the  bed,  can  be  obtained  by:  (a)  extrapolating, 
graphically  or  numerically,  the  £  vs.  G  curve  to  e  =  1.0  or  the  h  /h 
vs.  G  curve  to  h  /h  =  1.0  (Fan  and  Wen,  1961);  (b)  approximating  G   , 

paO  ^srr  ojjjgj. 

by  the  free  fall  terminal  velocity  of  a  single  particle  (Fan  and  Wen,  1961); 

(c)  using  the  following  semiempirical  equation  (Poddar  and  Dutt,  1969); 

1.687 

18  Re   .  +  2.7  Re   .     =  Ga  ,  (II-ll) 

ms I  3S  r 

where  Re   ^  is  the  Reynolds  number  at  the  maximum  semif  luidization  velocity; 
msf  '  ' 

(d)  using  any  of  the  empirical  correlations  by  Roy  and  Sarma  (1974a) ,  Roy 
and  Sen  Gupta  (1975)  or  Roy  (1971)  listed  in  Table  II-3. 
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PRESSURE  DROP 

The  overall  pressure  drop  of  both  fully  fluldlzed  and  semlf luidized 

beds  with  the  parameter  of  the  height  of  initial  static  bed  at  least  dense 

condition,  h  ,  has  been  experimentally  studied  by  Fan  and  Wen  (1961)  as 
o 

shown  in  Figure  II-4. 

Because  of  the  formation  of  the  packed  section  when  the  fluidized 
suspension  is  suppressed,  the  total  pressure  drop  should  be  the  algebraic 
sum  of  the  pressure  drops  across  the  fluidized  section  and  the  packed 
section,  since  these  two  sections  are  aligned  in  series  in  the  direction 
of  fluid  flow.  Therefore, 

AP  -  (4^),(h  ^  h   )  +  (4^)    h  (11-12) 

c     L  f      pa      L  pa   pa 

By  substituting  the  effective  weight  of  the  particles  into  the  first  term 
in  the  right-hand  side  of  Eq.  11-12,  the  Ergun  equation  (Ergun,  1952)  into 
the  second  term,  and  combining  the  resultant  expression  with  Eq.  II-l,  Fan 
and  Wen  (1961)  have  obtained 

(1  -e  )(hj  -h) 


i ]  a  -  Ej)(P3  -  pp 
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where  p,  and  fi   are  the  density  and  the  viscosity  of  the  fluid,  respectively, 

d  ,  the  equivalent  particle  diameter,  (p    ,  the  sphericity  of  particles,  u, 
p  s 

the  superficial  velocity  of  the  fluid,  and  G,  the  mass  flow  rate  of  the 
fluid  per  unit  area  of  the  empty  bed.   Equation  11-13  was  found  to  give  a 
reasonably  good  agreement  with  the  experimental  results  for  solid-liquid 
solid-gas  systems  with  particles  of  various  sizes  and  shapes  (Fan  and  Wen, 


1961;  Wen  £C  al. ,  1963).   However,  due  to  difficulties  in  determining  £^ 

and  e   accurately,  Eq.  11-13  often  led  to  values  below  those  measured 
pa 

(Wen  et  al. ,  1963;  Roy  and  Sen  Gupta,  1973;  Roy  and  Sarma,  1974b). 
Subsequently,  the  following  correction  factor  was  introduced  (Roy  and 
Sen  Gupta,  1973;  Roy  and  Sarma,  1974b); 

^ (11-14) 


^^^^Eq.  11-13 

This  correction  factor  was  empirically  correlated  by  Roy  and  Sen  Gupta 

(1973)  and  Roy  and  Sarma  (1974b)  as  a  function  of  the  system  variables, 

i.e. , 

D    p    h      h 

pre      o 

in  which  D  is  the  column  diameter,  h  ,  the  height  of  the  static  bed  under 
c  o 

the  least  dense  condition,  and  R,  the  bed  expansion  ratio.   These  corre- 
lations are  tabulated  in  Table  II -4. 

Similar  to  Eq.  11-13,  it  has  been  postulated  that  in  addition  to  the 
pressure  drop  for  the  fluid  flowing  through  packed  section,  the  fluid 
need  support  the  solid  particles  in  the  entire *bed  (Yutani  and  Fan,  1984). 
This  implies  that 


150(1  -  £   )^  (1  -  E   )   ^ 

(1  -  ep    1 

•  [(hf  -  h) ; r-]  —  (11-16) 

"  '■Ef   Epa-'   ^c 

It  appears  that  Eq.  11-13  tends  to  underestimate  the  pressure  drop  while 
Eq.  11-16  tends  to  overestimate  It;  the  actual  pressure  drop 
is  bounded  by  the  values  estimated  by  these  two  equations.  Therefore,  a 
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correction  factor,  c  ,  which  has  the  value  between  0  and  1,  may  need  be 
introduced  to  modify  them  as  (Yutani  and  Fan,  1984) 


AP^  =  (1  -  cp  hj(l  -  ep(p^  -  pp 

(1  -  e„  )(h  -  h) 

"  "^  "f '      c^f  -  4^      '  a-^fX'-s-V 
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MASS  TRANSFER 

The  mass-transfer  property  of  the  semif luidized  bed  was  Investigated 
by  Fan  et  al.  (1959;  1960),  using  the  benzoic  acid-water  system.   The 
performance  of  the  semifluidized  bed  as  a  mass-transfer  contactor  has  been 
found  to  be  a  linear  combination  of  those  of  the  packed  and  fluidized 
sections  in  the  bed.  The  experimentally  determined  overall  log-mean  mass- 
transfer  coefficient  can  be  expressed  by  the  correlation 

k.   -  CI  -  •'<)(k  )    +  X(k   )  (11-18) 

■Cm  r  LIB      pa  xm 

in  which  (k,)   and  (k.)   are  the  log-mean  mass-transfer  coefficients 
f  J,m       f  J,iii 

determined  for  the  fluidized  bed  and  packed  bed,  respectively.   Similar 
results  were  also  obtained  by  Govindrajan  et^   al.  (1976). 

ION  EXCHANGE 

Currently,  semifluidized  beds  have  been  widely  used  in  the  field  of 
ion  exchange.  Mobay  Chemical  Company,  Asahi  Chemical  Company  of  Japan, 
Bayer  AG  (Germany)  and  its  U.S.  counter  part,  Lewatit  Company,  are  advo- 
cating the  use  of  semifluidized  beds  with  ion  exchange  resins.   This  action 
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precipitated  from  the  discovery  that  a  fluidized  bed,  followed  by  a  fixed 
bed,  increased  the  efficiency  of  resin  utilization  by  improving  liquid- 
resin  contact.  The  fixed  bed  acts  as  a  polishing  section,  handles  the  ion 
leakage  from  the  fluidized  bed,  and  prevents  elutriation  of  resin  particles 
(Chementator,  1962;  1972;  Lewatit,  1972).   In  addition  to  a  higher  resin 
utilization  efficiency,  semifluidization  also  minimizes  the  voliraie  of 
regenerant  and  wash  water  needed,  reduces  pressure  drop,  and  operates  more 
consistently  than  the  conventional  fixed  bed  process. 

ADSORPTION 

Industrial  applications  of  adsorption  are  conducted  principally  in 
columns  of  fixed,  fluidized  or  moving  sorbent  beds.  Fixed  bed  adsorbers 
are  used  in  the  majority  of  industrial  applications  due  to  the  ease  and 
reliability  of  operation  of  these  units.  Mass  transfer  rates  are  relatively 
higher  in  fixed  beds,  but  pressure  drop  increases  markedly  as  the  particle 
size  is  decreased.  Fluidized  bed  operation  allows  the  use  of  relatively 
saaller  particles  and  provides  a  lower  pressure  drop  compared  to  packed 
beds. 

A  semif luidized  bed  adsorber  employed  to  remove  trace  organics  from 
aqueous  solutions  has  been  proposed  by  Mathews  and  Fan  (1983a;  1983b),  The 
adsorption  column  is  operated  as  a  semifluidized  bed  by  fluidizing  the  bed 
initially  and  then  compressing  the  fluid  bed  using  a  movable  sieve  plate. 
By  positioning  the  sieve  plate  at  various  heights,  different  ratios  of 
fluidized  bed  to  packed  bed  heights  can  be. obtained. 

Mathews  and  Fan  (1983a;  1983b)  compared  the  performance  of  fluidized, 
semifluidized  and  packed  beds  for  the  adsorption  of  phenol  from  aqueous 
solutions  onto  activated  carbon.   It  was  reported  that  the  same  amount  of 
adsorption  can  be  obtained  in  semifluidized  beds  as  in  packed  beds,  but  at 
a  lower  overall  column  pressure  drop,  Desorption  was  also  found  to  be 
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achievable  at  essentially  the  same  rate  as  in  packed  beds,  again  at  a 
lower  pressure  drop.   These  studies  indicate  that  solvent  recovery 
operations  involving  adsorption-desorption  operations  can  be  more  economi- 
cally achieved  in  semifluidized  beds.   In  their  earlier  work  (Mathews  and 
Fan,  1983a),  adsorption  studies  were  conducted  using  particles  of  a  constant 
size  at  several  fluidized  to  packed  bed  height  ratios  and  flow  rates,  and 
those  using  a  mixture  of  adsorbent  particle  sizes  were  also  conducted  by 
them  in  the  later  work  (Mathews  and  Fan,  1983b) . 

The  experimental  breakthrough  profiles  for  packed,  semifluidized  and 
fluidized  beds  and  the  pressure  drops  for  packed  and  semifluidized  beds 
presented  by  Mathews  and  Fan  (1983a)  are  shown  in  Figures  II-5  and  II-6, 
respectively. 

REACTOR 

Some  chemical  or  biochemical  reactions  can  be  carried  out  in  the 
semifluidized  bed  with  the  definite  advantage  over  the  operation  in  the 
fluidized  or  fixed  bed  reactor. 

Catalytic  Chemical  Reactor 

The  oxidation  of  benzene  has  been  studied  in  the  adiabatic  semlflui-: 
dized,  fluidized  and  fixed  bed  reactors  (Babu  Rao  and  Doraiswamy,  1970). 
The  catalyst  used  consisted  of  vanadium  and  molybdenum  oxides  on  silica  gel. 
Experiments  were  carried  out  at  an  inlet  temperature  of  380  °C  and  at  several 
benzene-air  ratios.  A  definite  improvement  has  been  observed  when  the 
oxidation  was  carried  out  in  the  semifluidized  bed  reactor  by  a  comparison 
of  the  conversions  obtained  in  three  types  of  reactors. 

Fixed  Film  Bioreactor 

The  use  of  the  semifluidized  bed  as  a  bioreactor  has  been  patented  by 
Fan  and  Wen  (1981;  1982).  The  idea  was  prompted  by  the  development  of  fixed 
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film  fluldized  bed  and  fixed  bed  bioreactors  (see,  e.g.,  Besik,  1973; 
Jeris  and  Owens,  1975;  Jewell  et  al. ,  1976;  1978;  Scott  and  Rancher, 
1976;  Weber  et  al. ,  1970).   In  the  fixed  film  fluidized  bed  bioreactor, 
the  fluidizing  particles,  usually  sand  or  activated  carbon,  are  coated 
with  viable  microbial  cells.  This  adhesion  of  microbes  onto  the  enormous 
surface  area  provided  by  the  fluidizing  particles  accounts  for  the  high 
concentration  of  biomass ,  Jeris  et  al.  (1977)  reported  the  concentration 
of  the  mixed  liquor  volatile  suspended  solids  (MLVSS)  in  the  fluidized  bed  - 
bioreactor  to  be  between  8,000  and  40,000  mg/jt ,  an  order  of  magnitude  higher 
than  that  in  a  conventional  suspended  growth  system.  Thus,  the  residence 
time  required  for  an  equivalent  task  should  be  one  to  two  orders  of  magnitude 
shorter  in  fixed  film  bioreactors  than  in  conventional  system  (Besik,  1973; 
Jeris  and  Owens,  1975;  Jewell  et  al . ,  1976;  1978).  The  use  of  semifluidized 
beds  in  place  of  fixed  and  fluidized  beds  in  biological  waste  treatment 
should  offer  distinct  advantages.  Among  these  are  the  elimination  of 
unstable  bed  expansion  and  elutriation  of  microbe  coated  particles  in  flui- 
dized bed  bioreactors,  and  the  reduction  of  bed  plugging  by  suspended  solids 
experienced  in  fixed  bed  operations.   Other  possible  advantages  of  using  a 
semifluidized  bed  bioreactor  in  place  of  a  fluidized  bed  bioreactor  have 
been  discussed  by  Fan  and  Hsu  (1980). 

The  patent  of  the  semifluidized  bed  bioreactor ,  proposed  by  Fan  and 
Wen  (1981) ,  describes  a  method  of  carrying  out  a  fluidized  bed  biological 
reaction  for  advanced  treatment  of  wastewater  in  which  the  top  of  the  bed 
is  restrained  with  a  media-retaining  perforated  plate  to  provide  a  smaller 
upper  fixed  bed  portion  and  a  larger  fluidized  bed  portion  therebeneath. 
The  wastewater  is  subjected  to  initial  treatment  in  the  fluidized  portion 
in  the  bed  and  to  final  treatment  in  the  fixed  portion  of  the  bed.  The 
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media-retaining  perforated  plate  together  with  the  filtering  action  of 
the  fixed  bed  portion  substantially  prevents  media  loss  and  greatly  reduces 
the  bacteria  loss.  In  addition,  more  complete  removal  of  the  organic 
contaminants  is  obtained  without  sacrificing  the  high-rate  capacity  of 
the  fluidiaed  bed. 

The  fixed  bed  section  drives  the  reaction  to  completion  because  of 
the  higher  blomass  concentration  thereon.  Viable  bacteria  cells  escaping 
from  the  fluidized  bed  are  trapped  by  the  filtering  and  retention  action 
of  the  fixed  bed  section.  Further,  the  entering  wastewater,  which  is  high 
in  contaminants,  is  first  treated  in  a  fluidized  bed  where  intensive  mixing 
occurs  while  the  effluent  from  the  fluidized  portion,  which  is  lean  in 
contaminants,  is  next  treated  in  a  fixed  bed  where  mixing  is  not  Intensive. 
This  arrangement  provides  more  efficient  waste  removal  for  given  volume. 

The  method  of  this  invention  also  facilitates  cleaning  of  the  reactor 
beds.   Periodically,  the  treatment  of  the  wasterwater  is  discontinued  and 
the  perforated  plate  is  raised  to  a  level  out  of  contact  with  the  bed  when 
fully  fluidized  and  above  an  outlet  which  is  normally  above  the  plate. 
Wash  water  is  then  passed  upwardly  through  the  bed  at  a  velocity  fully 
fluidizing  the  bed.  The  wash  water  is  removed  together  with  excess  cellular 
material  through  the  outlet,  which  is  preferably  provided  in  the  side  of 
the  reactor.  The  fixed  portion  of  the  bed  is  broken  up  to  release  the 
cellular  material  including  the  bacterial  floe.  This  may  be  promoted  by 
the  use  of  a  mechanical  agitator  within  within  the  fixed  bed  portion,  such 
as  a  shearing  mixer. 

Another  patent  of  the  semif luidized  bed  bioreactor,  also  proposed  by 
Fan  and  Wen  (1982),  suggests  a  new  method  for  cleaning  of  the  reactor  beds. 
During  the  cleaning  period,  the  perforated  plate  is  raised  to  a  level  at 
which  the  bed  can  be  fully  fluidized  but  below  a  wash  water  outlet  which 
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is  normally  above  the  plate.  Wash  water  Is  then  passed  upwardly  through 
the  bed  at  a  velocity  fully  fluidizing  the  bed.   The  wash  water  is  removed 
together  with  excess  cellular  material  through  the  outlet,  which  is  preferably 
provided  in  the  side  of  the  reactor.  The  fixed  portion  of  the  bed  is  broken 
up  to  release  the  cellular  material  including  the  bacterial  floe.   This  may 
be  promoted  by  the  use  of  concurrent  gas  introduction,  such  as  combined 
flow  of  water  and  air  under  high  turbulence. 

POROUS  MEDIA  FILTER 

The  use  of  the  semifluidized  bed  as  a  porous  media  filter  was  first 
patented  by  Wen  and  Fan  (1979) .   The  patent  describes  a  method  of  operating 
a  semifluidized  bed  filter  in  which  the  filter  cake  thickness  is  controlled 
by  the  fluidizing  action  of  the  porous  media  or  grains  instead  of  complicated 
mechanical  devices  such  as  cutting  blades, 

A  conventional  deep  bed  filter  is  composed  of  a  fixed  or  packed  bed 
of  grains,  usually  sand  or  sand  and  anthracite.  The  suspended  solids  in  the 
effluent  are  retained  in  the  filter  either  by  adhesion  to  the  filter  grains 
or  by  the  straining  process.   The  operation  of  a  deep  bed  filter  is  typically 
characterized  by  shallow  penetration,  rendering  much  of  the  bed  ineffective 
as  a  filtering  medium.   Furthermore,  a  layer  of  filtered  solids  will  collect 
on  the  surface  and  within  the  first  few  layers  of  the  grains,  giving  rise 
to  a  sharp  increase  in  the  pressure  drop. 

In  a  semifluidized-bed  filter  proposed  by  Wen  and  Fan  (1979)  ,  the 
packed  section  of  the  bed  functions  as  a  deep  bed  filter  whereas  the  fluidized 
section  provides  the  turbulent  scouring  action  in  delaying  clogging  of  the 
bed.  The  use  of  this  apparatus  for  filtration  purposes  presented  by  them 
is  shown  in  Figure  II-7.   Initially,  the  bed  will  be  in  a  completely  fixed 
configuration  as  the  fluid  stream  containing  suspended  solids  enters  the  bed. 
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The  suspended  solids  upon  entering  the  portion  of  the  bed  adjacent  to  the 
distributor  zone  will  be  unable  to  penetrate  the  bed  to  any  great  extent 
before  being  stopped  by  the  grains,  A  "cake"  of  suspended  solids  will 
begin  to  be  built  up  in  this  zone.  The  formation  of  the  "cake"  will  produce 
an  abrupt  increase  in  the  pressure  drop.   This  rise  in  the  pressure  drop 
will  activate  a  driving  mechanism  which  will  raise  the  upper  sieve  plate. 
The  mechanical  movement  of  the  sieve  plate  in  conjunction  with  the  forces 
being  exerted  by  the  fluid  stream  will  cause  disruption  of  the  formed  "cake" 
of  suspended  solids  into  small  fragments.  Cake  disruption  will  also  be 
aided  by  the  fluidized  particles.  Also  the  act  of  disruption  of  existent 
"cake"  will  cause  more  of  the  particles  which  are  encased  by  the  cake  to 
become  fluidized.   Raising  of  the  sieve  plate  will  continue  while  small 
suspended  solids  in  the  stream  are  being  trapped  in  the  fixed  bed.  The 
amount  of  fixed  bed  will  gradually  decrease  with  successive  disruption  of 
the  formed  cakes.  The  filtration  cycle  is  terminated  when  the  suspended 
solids  penetrate  through  the  ever  decreasing  packed  bed  section. 

A  recent  study  of  semifluidized-bed  filters  was  presented  by  Hsu  and 
Fan  (1981;  1982),  and  Fan  and  Hsu  (1984).   In  their  work,  the  filtration 
cycle  starts  with  the  restraining  plate  in  the  proper  position  forming  a 
semifluidized  bed  with  25%,  50%  or  757=  fluidization.   Instead,  of  forming 
the  "cake"  around  the  entrance  section  of  the  fixed  bed,  the  suspended 
solids  together  with  some  of  the  fluidizing  grains  will  begin  to  form  a 
fixed  bed  in  one  side  of  the  initially  fluidized  section,  thus  lowering 
the  sharp  increase  in  the  pressure  drop.  With  the  filtration  run,  the 
newly-built  fixed  bed  will  extend  gradually  down  to  the  distributor  plate. 
Their  results  indicate  that  the  filtration  run  duration  can  be  increased  up 
to  6  times  that  of  a  deep  bed  filter  with  the  use  of  a  semifluidized  bed 
without  sacrificing  the  filtrate  quality. 
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HEAT  TRANSFER 

Purachander  Rao  and  Kapathi  (1969)  studied  the  variation  in  the 
magnitude  of  wall  heat-transfer  coefficient  with  changes  in  the  operating 
parameters  in  a  semif luidized  bed.  Their  experimental  results  show  the 
values  of  the  Nusselt  number  obtained  for  the  semifluidized  bed  to  be 
lower  than  those  for  the  fixed  bed  and  higher  than  those  for  the  fluidized 
bed  at  a  given  Reynolds  number.   It  was  also  shown  that  the  values  of  the 
heat-transfer  coefficient  for  the  semifluidized  bed  can  be  calculated  from 
the  available  correlations  for  the  heat-transfer  coefficients  in  packed 
beds  and  fluidized  beds  (see,  e.g..  Leva  et   al. ,  1948);  the  mean  heat- 
transfer  coefficient  of  the  semifluidized  bed  was  shown  to  be  the  sum  of 
the  heat-transfer  coefficient  of  the  packed  and  fluidized  beds  multiplied 
by  the  length  fraction  of  the  respective  sections, 

SEMIFLUXDIZATION  WITHOUT  RESTRAINING  PLATE 

The  upper  restraining  plate  is  not  essential  to  the  formation  of  a 
semifluidized  bed.  A  new  type  of  solid-gas  semifluidized  bed  contactor 
was  proposed  by  Corella  and  Bilbao  (1982).  The  contactor  was  constructed 
from  a  column  with  a  considerable  enlargement  in  the  cross  section  of  its 
upper  zone.  This  enlargement  in  cross  section  caused  a  reduction  in  the 
gas  velocity  in  that  region.  The  inlet  gas  was  introduced  in  such  a  way 
that  the  gas  velocity  was  greater  than  the  minimum  fluidization  velocity 
in  the  lower  zone  and  lower  than  the  minimum  fluidization  velocity  in  the 
.upper  zone,  thus  giving  rise  to  the  formation  of  a  fluidized  bed  and  a  fixed 
bed  in  series.  The  fluid  dynamics  and  bubble  behavior  in  this  type  of 
contactor  was  also  investigated, 

SEMIFLUIDIZATION  IN  A  GAS-LIQUID-SOLID  SYSTEM 

Fluidization  in  gas-liquid-solid  systems  where  liquid  is  the  continuous 
phase  can  be  achieved  with  either  cocurrent  upflow  of  gas  and  liquid  or 
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downflow  of  liquid  counter  to  the  upflow  of  gas.  The  downf lowing,  liquid- 
supported  bed  is  termed  an  inverse  gas-liquid-solid  fluidized  bed.  A  semi- 
fluidized  bed  is  formed  when  a  mass  of  fluidized  particles  is  compressed 
from  above  with  a  porous  restraining  plate. 

The  hydrodynamics  of  inverse  and  cocurrent  gas-liquid-solid  semi- 
fluidization  have  been  experimentally  investigated  by  Chem  et  al.  (1983; 
1984).   Their  investigation  was  concerned  with  the  behavior  of  the  pressure 
drop,  gas  holdup,  onset  liquid  velocity  for  semifluidization,  and  heights 
of  the  packed  section  and  fluidized  section. 
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NOMENCLATURE 

A  "  cross-sectional  area  of  the  bed 

C  =  correction  factor  defined  in  Eq.  11-14 

C-  =  correction  factor  used  in  Eq.  11-17     ;        ^  ^ 

d  =   equivalent  particle  diameter 
P 

D  =  column  diameter 
c 

g  =  conversion  factor 

G  =  mass  flow  rate  of  fluid  per  unit  area  of  the  empty  bed 

G  £  =  onset  velocity  of  fluidization 

G   ,  =  onset  velocity  of  semif luidization 
osf  -^ 

G  £  =  maximum  semif luidization  velocity 
Ga 


3  2 

^   Pf(p  -  Pf)g/)Lt  J  Galileo  number 


P 
h  =  semif luidized  bed  height 

fe  =  height  of  initial  static  bed  at  the  least  dense  condition 

0 

h-         =  height  of  fully  fluidized  bed 

h  =  height  of  packed  section  in  semif luidized  bed 

k*         =  overall  log  mean  mass-transfer  coefficient 

(k_)  „       =  log  mean  mass-transfer  coefficient  of  the  fluidized  section 
in  the  semif luidized  bed 

(k  )p      =  log  mean  mass-transfer  coefficient  of  the  packed  section  in 
the  semif luidized  bed 

AP.         =  overall  pressure  drop  through  semif luidized  bed 

^')f  ~   pressure  drop  per  unit  length  in  the  fluidized  section  of  a 

semif luidized  bed 
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AP 
(— )         =  pressure  drop  per  unit  length  in  the  packed  section  of  a 

semifluidized  bed 

R  =  bed  expansion  ratio,  semifluidized  bed  height  divided  by- 

initial  static  bed  height  under  the  least  dense  condition 

RiB   .        =  (d  G  ^)/M»  Reynolds  number  at  the  maximum  semif luidization 
SOSC  p  msf 

velocity 

Ra   -         =  (d  G   ^) /u,  Reynolds  number  at  the  onset  of  semif luidization 
osf  p  osf  ^ 

3 
Sf  =(W-W  )/[(h-h)p],  semif luidization  number 

u  =  superficial  fluid  velocity 

u  =  terminal  velocity  of  the  particle 

U  -  =  minimum  fluidization  velocity 

ffix 

W  =  total  weight  of  particles 

W  =  weight  of  particles  in  the  packed  section  of  a  semifluidized 

P*  bed 

X  =  weight  fraction  of  particles  in  the  packed  section  of  a 

semifluidized  bed 

g  =  porosity  of  the  initial  static  bed  under  the  least  dense 

dense  condition 

e-  =  porosity  of  the  fluidtzed  section  or  the  fully  fluidized 

bed 

£  =  porosity  of  the  packed  section 

^  =  sphericity  of  particles 

fl  =   viscosity  of  fluid 
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Figure  II-2.   A  semif luidized  bed  in  which  the  density  of 
particles  is  less  than  that  of  the  fluid. 
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Figure  II-3.  Comparison  of  measured  heights  of  packed  section 
with  theoretical  values  for  50x60  mesh  and  40x45 
mesh  particles  (Fan  and  Wen,  1961). 
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Figure   II-4.      Overall   pressure  drop   through  fully  and   semif luidized 
beds   of  50x60  mesh  particles    (Fan  and  Wen,    1961). 
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Hydraulic   Loading  Rate, m/min 


Figure  II-6.  Pressure  drop  data  for  F-400  carbon: 

(O)  packed  bed  and  semif luidized  beds  with 
(D)  25%  fluidization,  (•)  50%  fluidization  and 
(■)  75%  fluidization  (Mathews  and  Fan,  1983a). 
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(a)  Initial  Stage. 


(b)  Cake  of  Fines  Start  to  Form  at 
the  Bottom  of  the  Fixed  Bed. 


r 


M\ 


m- 


(c)  Top  Sieve  Plate  is  Gradually  (d)  Final  Stage  of  Filtration. 
Lifted  and  a  Fluidized  Bed  is       Fines  Collected  Can  Be  Purged 
Formed  Making  Flow  and  Filtra-      by  Lifting  the  Sieve  Plate  and 
tion  Easier.  Washed  by  Suitable  Fluid. 


Figure  II-7.  Use  of  the  Semif luidized  Bed  for  the  Filtration  Process 
(Wen  and  Fan,  1979) 


CHAPTER  III  mSTER  EQUATION  AKD  ONE-STEP  PROCESSES 

Macroscopic  physics  deals  with  macroscopic  quantities,  such  as  the 
positions  and  velocities  of  bodies,  electrical  charges  and  currents, 
amounts  or  concentrations  of  chemical  compounds.  They  obey  macroscopic 
laws  which  can  be  derived  from  general  principles  such  as  conservation 
laws,  together  with  some  specific  assumptions  of  phenomenologtcal  nature, 
for  example,  those  of  Fourier  and  Fick.  The  resultant  differential  equa- 
tions are  the  macroscopic  or  deterministic  equations;  for  instance,  the 
rate  equations  for  chemical  reactions,  the  damped  harmonic  oscillator 
equation  for  a  pendulum  in  air,  the  diffusion  equation.   They  are  approxi- 
mate; in  reality  there  are  small  deviations  from  them,  which  show  up  as 
fluctuations . 

Of  course,  the  fluctuations  cannot  be  found  exactly;  that  would  be 
tantamount  to  solving  the  microscopic  equations.   In  principle,  all  informa- 
tion is  contained  in  the  microscopic  equations  of  motion  of  all  particles, 
but  It  hardly  needs  saying  that  an  exact  solution  of  these  equations  Is 
beyond  human  means.   It  is  therefore  necessary  to  embark  upon  a  less 
ambitious  program  and  develop  a  theory  which  goes  beyond  the  macroscopic 
description  in  that  it  Includes  fluctuations,  but  shortcuts  the  connection 
with  the  microscopic  equations.  This  is  the  so-called  stochastic  approach 
which  views  the  macroscopic  phenomena  as  the  outcome  of  the  collective 
behavior  of  numerous  particles  and  includes  fluctuations  by  means  of  the 
repeated  randomness  assumption;  some  call   It  the  mesoscopic  level  of 
approach. 

One  popular  stochastic  approach  starts  out  from  the  transition  proba- 
bility density,  P^,  and  assxjmes  that  it  obeys  the  Chapman-Kolmogorov  equation 
with  the  type  of 
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^  P  =  P  U  (III-l) 

at  —   —  — 

where  U  is  a  matrix  of  intensities  of  transition, 

MASTER  EQUATION 

The  Master  equation  is  an  equivalent  form  of  the  Chapman-Kolmogorov 
equation  for  Markov  processes,  but  it  is  easier  to  manipulate  and  more 
directly  related  to  physical  process  to  be  modeled. 

Consider  a  discrete  set  of  states  with  labels,  n,  the  Master  equation 
can  be  written  as  (see,  e.g.,  Van  Kampen,  1981) 


=  2.  [V  ,„  P  .(t)  -  v„  ,  p  (t)]  (III.2) 


this  form  the  meaning  becomes  particularly  clear.   The  Master  equation 
is  a  gain-loss  balance  equation  for  the  probability  of  being  in  state  n; 
the  first  term  in  the  right-hand  side  of  this  expression  is  the  gain  due 
to  transitions  from  other  states,  and  the  second  term  is  the  loss  due  to 
transitions  into  other  states. 

As  an  example  for  demonstrating  the  usefulness  of  the  Master  equation, 
let's  consider  a  simple  one-step  or  birth-death  process.   For  the  birth- 
death  process,  it  is  assun 
(see,  e.g.,  Chiang,  1980) 

1.  The  conditional  probability  that  a  birth  event,  which  changes 
the  state  from  n  to  (n+1) ,  will  occur  during  the  time  Interval 
(t,  t+At)  ±3    X      t  +  o(At),  where  X   is  the  intensity  of  birth 
transition. 

2.  The  conditional  probability  that  a  death  event,  which  changes 
the  state  from  n  to  (n-1) ,  will  occur  during  the  time  interval 
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(t,  t+At)  is  jx     t  +   o(At),  where  /j.     is  the  intensity  of  death 
transition. 
3.  The  conditional  probability  that  more  than  one  event  will  occur 
in  this  time  interval  is  o(iit),  where  o(At)  signifies  that 

lim  £^  =  0 

At>0   "^ 

Substituting  the  corresponding  )(,  's  and  n   's  into  the  Master  equation,  Eq. 
III-2,  and  neglecting  the  terms  with  the  order  of  At,  i.e.,  o(At),  give 

%^  =  ^n-1  Vl(')  -^/^n+l  Vl('>  -  ^K^l'r?    ^n^'^  ' 

(III-3) 


-it-  =-^Po(t>+^i  Pi(^>        ■  ("^-^> 

aod 

dp  (t) 

-. V^V^(^)-/^n„Pn„<0  aiI-3) 

For  convenience,  define  ^»,  J  -i  »  ^   ^nd  ^  ^  to  be  zero.   Equations  III-3 
through  III-5  can  therefore  be  combined  as 

4i-  =  ^n-l  Vl^')  ^/^n+l  Pn+l^t)  "  ^^n  ^^^n^  Pn^'^ ' 

n=0,  1,  2,  ....  ng         (III-6) 

The  solution  of  this  expression  gives  the  probability  distribution  of  the 
process. 

ONE-STEP  PROCESSES  AM)  ONE-STEP  OPERATOR 

The  one-step  or  birth-death  processes  are  a  special  class  of  stochastic 
processes,  which  occur  in  many  applications  and  can  be  analyzed  in  some  detail. 
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The  one-step  operator,  E»  is  defined  by  its  effect  on  an  arbitrary 
function  f(n)  as  (see,  e.g.,  Van  Kampen,  1981) 

Ef(n)  =  f(n+l)  and  E"^f(n)  =  f(n-l)  (III-7) 

If  n  ranges  from  -oo  to  +«  ,  it  is  possible  to  regard  E  as  an  actual  operator 
In  function  space,  but  if  there  are  one  or  two  boundaries,  it  is  better  to 
regard  E  simply  as  an  abbreviating  symbol.   Most  of  its  properties  are 
immediately  evident,  but  we  note  in  particular 

£  g(n)Ef(n)=  2°  f(n)  E'-^g(n)  (III-8) 

n=0  n=l 

for  any  pair  of  functions  f  and  g.  The  difference  in  the  limits  of  the 
summations  in  both  members  is  inconvenient.   Fortunately,  it  is  often  useful 
in  practice,  either  because  the  terms  are  infinite  (and  f  or  g  vanish  suffi- 
ciently rapidly  at  infinity)  or  the  missing  terms  may  be  added  to  the  sum 
since  they  are  zero.   For  instance,  the  term  with  n=0  may  be  added  on  the 
right-hand  side  of  this  equation  if  f(0)  or  g(-l)  is  defined  to  be  zero. 
Therefore,  this  equation  becomes 

i  i" 

2°  g(n)  Ef(n)  =  2°  f(n)  E'^g(n)  (III-9) 

n=0  n=0 

Equation  III-9  means  that  f  and  g  in  the  summation  can  be  exchanged  by 
changing  the  sign  of  the  one-step  operator. 

With  the  aid  of  the  one-step  operator,  the  Master  equation  for  one-step 
processes,  Eq.  III-6,  may  be  rewritten  as 

dp^(t) 

-5S—  =  (E-  1)  tx^   p^(t)  +  (E   -  1)  \^   P^(t);  (III-IO) 
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With  the  aid  of  Eq.  III-9,  the  mean  and  variance  of  n  governed  by  this 

equation  can  be  easily  obtained.  As  an  example,  assume  that 
intensities  of  transition  take  the  form  as 

X^   =  a(ng  -  n)  (III-ll) 

ix  =  fln  (111-12) 

Substitution  of  Eqs.  III-ll  and  III-12  into  Eq.  III-IO  yields 

dp„(t)  1 

-jS—  =  (E  -  1)  /3n  p^(t)  +  (E   -  1)  a(ng  -  n)  p^(t)i   ; 

n  =  0,  1,  2,  ...,  ng       (III -13) 

Multiplying  both  sides  of  this  expression  by  the  respective  n's  and  summing 
all  resultant  equations ,  we  have 

■    ^=  2°  n(E  -  1)  ^n  p^(t) 
n=0 

+  2°  n(E'^  -  1)  0(^0  -  ")  PnC^)         (III-14) 
n=0 

where  <n>   is   the  mean  of  n.      From  the  property  of  the  one-step  operator  as 
given  in  Eq.    III-9,  we  see   that  Eq.    111-14  is  equivalent  to 


d<n>  —Q  ,   .  ,_-l        . , 

-rr  =       2    /3n  p„(t)(E       -  l)n 


+    2°  a(n.   -n)    p„(t)(E-   l)n  (III-15) 


=  -^n>  +  a(n     -  <n>) 

The  solution  of  this  equation  subject  to  the  initial  condition, 

<n>  =0   at  t  =  0, 

gives 

<n>  =^[1  -e-("+«t  (III-16) 
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The  variance  of  n,    <<n>>,  defined  as 


is  also  of  interest.      Multiplying  both  sides   of  Eq.    III-13  by  the   respective 


^=       2°n'(E-  l)^np^(t) 
n=0 


+    2°  n^E"^   -   1)   aCng   -  n)    p^(t)  (III-17) 

n=0 


From  Eq.  III-9,  the  above  equation  reduces  to 


d<n  >  _   0      /^\/i."^   1%  2 
-J^-     2  ^nP„(t)(E   -  l)n 
n=0 


+  2°  a(nr,  -  n)  p„(t)(E-  l)n^ 


=  -2(a+^)<n^>  +  [(2ng  -  l)a  +/3]<n>  +  an^  (IIX-18) 

The  solution  of  this  expression  subject  to  the  initial  condition, 

2 
<n  >  "  0  at  t  =  0, 

yields  the  variance, <<n>>,  as 

«"»  =  5+S  [^  -  ^    "']{!-  '■    „^o '■}  (III-19) 
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OTHER  APPLICATION  OF  THE  ONE-STEP  OPERATOR  TECHNIQUE 

The  one-step  operator  technique  is  not  restricted  to  solve  the 
linear,  one-step  processes;  it  can  be  used  elsewhere.  As  an  example 
consider  a  simple  polymerization  system  (Fan,  1983) .   Let  M  be  the  monomer, 
P.  be  the  free  radical  with  j  chain  length,  P,  be  the  polymer  with  j  chain 
length.   The  polymerization  mechanism  is  assumed  to  be 
k. 


M  ^P^ 

;  initiation 

k 

p-  +M— E-* 

k 
P-  +  M— 2^ 

'J+1 

;  propagation 

Vi 

;  termination 

The  corresponding  batch  kinetic  equations  are 
^=-k^[M]  -  (k+kJ[M][P^] 


(III-20) 


k. [M]    -    (k  +k  )[M][P:] 
1  p     m  i 


d[p; 


-^^  =  kp[M][p-.^]    -    (kp+kj[p-]      ;    j>2 


d[P.l 
dt 


k  [Mi[p;    ] 


;  2^2 


(III-21) 


(III-22) 


(III-23) 


[P,]   =  ^i,    [Pj] 


By  means   of   Che  eigenzeit  transformation  which  defines 
dT  =    [M]    dt    , 
Eqa.   III-20  through  III-23   can  be   transformed,    respectively,    as 


d[M]  ,       , 


(III-24) 


(III-25) 


d  [p  ■  1 

-^-^  =  k.    -    (k  +k  )  [p:  ]  (III-26) 

dT  1        ^   p     m       1 

d[p-.] 

-T-J-  =  k^[p;   J    -    (k  +k  )[p:]    ;    j>2  (III-27) 

dt  p      J-I  p     m        J 


d[P.] 

-T-J-  =  k    [p;      ]  ;    ji2  (III-28) 

dT  m     j-i 

Notice  that  these  are  the  linear  equations. 

Addition  of  Eqs .  III-26  and  III-27  leads  to 

d[P:] 

-r-i-  =  k  -  k  [PM  (III-29) 

qT      1    m  T 

The  solution  of  this  equation  subject  to  the  initial  condition 

0, 


yields 


k       -k  T 
[P^]  =  j^  (1  -  e  "  )  (III-30) 


The  mean  and  variance  of  the  chain  length  of  polymer  can  be  related 
to  the  first  and  second  moments  of  the  MWD  (molecular  weight  distribution). 
The  general  expression  is 

CO 

E(j"]  =  .!,_  j"{[P-]  +  [P*]  }       -■■-    -   -  •  ■'     (III-31) 

2 
where  E[ j]  is  the  first  moment  and  E[ j  ]  is  the  second  moment.   The  technique 

used  in  the  preceding  section  can  also  be  employed  here  to  obtain  the 

moments  defined  by  Eq.  III-31.   Multiplying  both  sides  of  Eqs.  III-26  through 

III-28  by  respective  j 's  and  summing  all  resultant  equations  result  in 


dE[ 
dt 


Introducing  the  one-step  operator,  Ej  this  expression  reduces  to 


From  Eq.  III-9,  we  see  that  the  above  equation  is  equivalent  to 

dE[ 


(III-33) 


fJ-U-  =  k.  -  (k  +k  )  .2,  [p;  ](i  -  E)j 

dt      1   ^  p  m  j=l   j 

00 

k  +  (k+k  )   2  [p;] 

1     p  m  j-1    J 


=  k.  +  (k  +k  )  [P' 1  (III-34) 

i   ^  p  m'  ■■  T  ' 

Substituting  Eq.  III-30  into  Eq.  III-34,  the  solution  of  the  resultant 

equation  subject  to  the  initial  condition 

,   i 

E[j]  =  0  at  T  =  0,  .    ' 

gives 

k,  (k  +Zk  )  k,  (k  +k  )        -k   T  '        '        ' 

E[jl=    ['I       "    ]x+    [^      P     °    ]e     -^ 


k.  (k  +k  ) 
"     P,  "  (III-35) 


The  mean  is  equal  to  the  first  moment. 

The  derivation  of  the  second  moment  follows  the  similar  procedure. 

2, 
Multiplying  both  sides  of  Eqs.  III-26  through  III-28  by  respective  j   s 

and  summing  all  resultant  equations  give 
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dx  1        ^   p     m     j=i  J 


=  k.  -  (kp+kj  .|  j^i-e-^[p:] 

=  k.   +   (kp+kj{2   jlj^   J[P-]   +  [P^]}  (iri-36) 

00 

To  solve  the  above  equation,  .2  j[P*]  has  to  be  found.   Multiplying  both 
sides  of  Eqs.  III-26  and  1X1-27  by  respective  j 's  and  summing  all  the 
resultant  equations  result  in 

d     1^    j[p:]  a  CO 

■^"dt    ^  =  ^i  -  VV  j5i  j[^-]  +  s  jii  j[p-.i]       '   ~  ■ 

Substituting  the  solution  of  this  equation  and  Eq.  III-30  into  Eq.  III-36, 
the  solution  of  the  resultant  equation  subject  to  the  initial  condition 


yleM 


E[j^]  =0  atT=  0, 


,    2k,k  (k  +k  )   -k  T    k.  (k  +k  )  (4k  +3k  )   -k  t 
E[j']  =[   ^  P^  P  -^  ]xe  "   +[  ^  P  °3   P   ^  ]e  " 


k.  (k  +k  )  (2k  +3k  )  +  k  ^      k,  (k  +k  )  (4k  +3k  ) 

+[  "   P  "' E-^2 HL_],  .  ^  P  -    P   -    (III.38) 

k  ^  k  "^ 


The  variance  can  be  calculated  from  the  first  and  second  moments. 


CHAPTER  IV  ANALYSIS  OF  DEEP  BED  FILTRATION  DATA: 
Modeling  as  a  Birth-Death  Process 

In  this  chapter,  a  fairly  general  stochastic  model,  namely,  the  linear, 
birth-death  process,  is  employed  in  conjunction  with  the  Carman-Kozeny 
equation  to  simulate  the  performance  of  the  deep  bed  filter  in  terms  of  the 
pressure  drop  dynamics  under  a  constant  flow  condition.   The  model  takes 
into  account  both  blockage  of  pores  by  suspended  particles  and  scouring  of 
deposited  particles;  the  entire  bed  is  assumed  to  be  spatially  lumped.   In 
spite  of  this  simplicity,  the  model  appears  to  represent  the  majority  of  the 
available  experimental  data. 


A  stochastic  birth-death  model  was  first  employed  by  Lltwlniszyn  (1965) 
to  describe  pore  blockage  and  scouring  in  a  filtration  process.   By  incorpo- 
rating the  Carman-Kozeny  equation,  this  model  has  been  developed  to  simulate 
the  pressure  drop  dynamics  through  the  deep  bed  filter  under  the  constant 
flow  condition. 

Birth-Death  Process 

A  stochastic  process  {N(t)}  is  a  family  of  random  variables  describing 
an  empirical  process  whose  development  is  governed  by  probabilistic  laws. 
In  considering  the  filtration  process  in  a  deep  packed  bed  as  a  stochastic 
process,  the  number  of  blocked  pores  in  a  unit  volume  of  the  bed  at  time  t, 
N(t) ,  can  be  taken  as  the  random  variable  (Litwiniszyn,  1966);  a  specific 
value  of  N(t)  will  be  denoted  by  n.  Given  N(t)  =  n  it  is  assumed  that  for 
the  birth-death  process  (Chiang,  1980),  ''  "...   . '"  ' 

1.  The  conditional  probability  that  a  birth  event  will  occur,  i.e. , 
that  an  open  pore  will  be  blocked  during  the  time  interval 
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(t,  t  +  At) ,  is  X  AC  +  o(At),  where  X  Is  a  function  of  n. 

2.  The  conditional  probability  that  a  death  event  will  occur,  i.e., 
that  a  blocked  pore  will  be  scoured  during  the  time  interval 

(t,  t  +  At),  is  u  At  +  o(At),  where  y^  is  a  function  of  n. 

3.  The  conditional  probability  that  more  than  one  event  will  occur  In 
this  time  Interval  is  o(At),  where  o(at)  signifies  that 

Obviously  the  probability  of  no  change  in  the  time  interval  (t,  t  +  At)  Is 

n  -  X  At  -  u  At  -  o(At)] .   The  probability  that  exactly  n  pores  are  blocked 
■^     n      n 

at  the  moment  t  will  be  denoted  as  P^(t)  =  Pr[N(t)  =n];  n=0,  1,  ... 
For  the  two  adjacent  intervals  (0,  t)  and  (t,  t+At) ,  the  occurrence  of 
exactly  n  pores  being  blocked  during  the  time  Interval  (0,  t+At)  can  be 
realized  in  the  following  mutually  exclusive  ways. 

1.  All  n  pores  are  blocked  in  (0,  t) ,  and  none  in  (t,  t+At),  with 
probability  p  (t) [1  -  X  At  -  y^At  +  o(At)]. 

2.  Exactly  (n-1)  pores  are  blocked  in  (0,  t) ,  and  one  pore  Is 
blocked  in  (t,  t+At),  with  probability  Pj,_i  (t)  [x^_j^AC  +  o(At)]. 

3.  Exactly  (n+1)  pores  are  blocked  in  (0,  t) ,  and  one  pore  Is 
scoured  in  (t,  t+At),  with  probability  P„+x  ( ')  [P^+x'^'^  "*"  °(At)]. 

4.  Exactly  (n-j)  pores  are  blocked  In  (0,  t) ,  and  j  pores  are 
blocked  in  (t,  t+At),  where  2<jin,  with  probability  o(At). 

5.  Exactly  (n+j )  pores  are  blocked  in  (0,  t) ,  and  j  pores  are 
scoured  In  (t,  t+At),  where  2 < J  <  (n  -n) ,  with  probability  o(At). 

Considering  all  these  possibilities  and  combining  all  quantities  of  o(At) 
give 
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+  P  ,,(t)M  _„Ai:  +  o(At),   n  >  1  (IV-1) 

n+i    n+i  _~ 

and 

Pg(t+AC)   =   pg(t)[l  -  AgAtl  +  p^(t)p^At  +  o(At)  av-2) 

Rearranging  these  equations  and  taking  the  limit  as  At  -*  0  yield  the 
following  Master  equations  (see  APPENDIX  A) 

dp  (t) 

—2 =  A   ,p  At)    -    i\     +  M„)P„(0  +  M  ,,P  ,,(t),   " 

dt         n-1  n-1        n    n   n       n+i  n+i 


(IV-3) 


and 


Litwiniszyn  (1966)  has  assumed  that  the  intensities  of  transition, 
A   and  U  ,  take  the  form 


X^   =   a(ng  -  n)  ,    n=  0,1,2,. ..,n^  '  (IV-5) 

U^   =   Bn  ■  ,  'r.  (IV-6) 

where  n  is  the  total  number  of  open  pores  susceptible  to  blockage  at  the 

0 
moment  t  =  0,  and  a  and  g  are  the  proportionality  constants.   The  constants 

a  and  g  may  be  called  the  blockage  and  scouring  constants,  respectively. 

Equation  IV-5  implies  that  the  rate  of  pore  blockage  is  proportional  to 

the  number  of  open  pores,  and  Eq.  IV-6  Implies  that  the  rate  of  pore 

scouring  is  proportional  to  the  number  of  blocked  pores.   Introducing  Eqs. 


IV -4 

IV-5  and  IV-6  into  Eqs.  IV-3  and  IV -4,  respectively,  the  following 
equations  are  obtained; 

dp  (t) 


dt 


aed 


dc 


-«ngPg(t)     +     3p^(t) 


(IV-7) 


(IV-8) 


At  the  start  of  the  filtration  process,  all  pores  are  open;  thus,  the 
initial  conditions  to  Eqs.  IV-7  and  IV-8  may  be  expressed  as 

P^(0)   =   0     n  =  1,2 ng  (IV-9) 


Solution  of  Eqs.  IV-7  and  IV-8,  subject  to  initial  conditions,  Eq . 
IV-9,  yields  the  distribution  of  the  probabilities.   However,  such  infor- 
mation has  limited  physical  significance.   Instead,  the  mean  number  of 
blocked  pores  at  a  given  moment  t,  namely,  E[N(t)],  is  determined;  it  is 
defined  as  ' 

"0 
E[N(t)]   =    j;   np  (t)  (IV-10) 

n=0   " 

To  evaluate  this  value,  it  is  unnecessary  to  solve  Eqs.  IV-7  and  IV-8. 
We  can  use  the  method  of  probability  generating  function  defined  as 

G(s,t)   =     1°   s"p  (t)  (IV-U) 


Multiplying  both  sides  of  Eq.  IV-7  by  the  respective  s  's  and  Eq.  IV-8 
by  s  ,  and  summing  all  resultant  equations,  the  following  expression  is 
obtained  upon  rearrangement  (see,  e.g.,  Chiang,  1980;  also  see  APPENDIX  B). 


1G(S^   ^   3G(s,t)^g  ^  (^  .  g)^  .  ^^2j 
3t         3s 


(IV-12) 


The  corresponding  initial  and  boundary  conditions  are 
G(s,0)   =   1 


(IV-13) 


G(l,t)   =   1 


(IV-14) 


The  solution  of  Eq.  IV-12  is 


G(s,c)   -   [ 


(as  +  6)  -  a(s  -  l)e 


-(a  +  6)t   no 


Examination  of  Eqs .  IV-10  and  IV-11  shows  that 

ElN(t))   -  i%^ 

s  =  1 

Substitution  of  Eq.  IV-15  into  Eq.  IV-16  gives 

-(a  +  B)t 


E(N(C))   =   on^I 


1  -  I 


0'     a  +  6 
The  variance  of  the  random  variable,  Var[N(t)],  takes  the  form  of 


(IV-15) 


(IV-16) 


(IV-17) 


C(s.t)   ^  8G(s,t)  _   3G(s,t)  2 


Var(N(t)]   =   (  "  -^-;"   +  ""'-'"^  -  [""^'■^M'-}| 

33  s  =  1 


7^11-  e-'-^^SJCj  ^^  .H^z_l 


■(a  +  B)c, 


'-}     (IV-18) 


Constant  Flow  Filtration 

To  utilize  the  results  obtained  in  the  preceding  subsection,  it  is 
necessary  to  relate  the  mean  of  the  number  of  blocked  pores  to  measurable 
variables.   In  the  present  modeling,  it  is  assumed  that  the  "pore"  refers 
to  finite  but  small  space  within  the  filter  bed,  which  is  susceptible  to 
the  deposition  of  solid  particles.   Each  "pore"  may  or  may  not  have  the 
clearly  defined  boundary,  and  the  space  surrounded  by  several  grains  or 
particles  may  be  composed  of  more  than  one.  pore.   Once  occupied  by  the 
particles  deposited  from  the  flow,  thus  preventing  additional  suspension 
to  pass  through  it,  the  pore  is  said  to  be  "blocked"  without  regard  to 
the  nature  of  deposition.   It  is  also  assumed  that  all  pores  have  the  same 
radius,  r.   For  the  constant  flow  filtration,  denoting  the  linear  velocity 
of  the  suspension  flowing  through  the  pores  at  the  moment  t  as  v(t)  and 
the  linear  velocity  at  the  onset  of  filtration  as  v  ,  we  have 

"O^^^O   °   dig  -  E[N(c)l)nr^v(t)  (IV-19) 

or,  upon  rearrangement, 

"''''>       '      Hq  -  E(N(t)l  "O  (IV-20) 

The  pressure  drop  for  laminar  flow  through  a  small  section  of  the 
packed  bed  can  be  adequately  described  by  the  Carman-Kozeny  equation  (see, 
e.g.,  Zenz  and  Othmer,  1950)  ' 

-  6P  =  150  i^^  ^  S.  (IV.21) 

£      d 

p 

The  superficial  velocity,  u,  in  this  equation  can  be  expressed  in  terms  of 


linear  velocity,  v,  as 
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Substitution  of  Eq.  IV-22  Into  Eq.  IV-21  yields 


(IV-22) 


(IV-23) 


P 
Thus,  the  pressure  drop  through  the  entire  packed  bed  becomes 

L     ,      , 
-  AP  =  /  150  (J^^^     ^^   dz  (IV-24) 

0       ^     d 
P 

For  simplicity,  the  entire  bed  Is  assumed  to  be  spatially  uniform 
and  single  set  of  values  for  e  and  v  is  employed  to  represent  the  pressure 
drop  characteristics  of  the  entire  bed.   Therefore,  Eq.  XV-24  reduces  to 

-f  =  150  (-4^)' -'"7  (IV-25) 

d 
P 

Note  that  the  change  of  fluid  velocity  in  the  pores  during  a  filtration 
run  is  expressed  by  Eq.  IV-20.   Thus,  substituting  Eq.  IV-20  into  Eq.  IV-25 
results  in  (see,  e.g.,  Hsu  and  Fan,  1984). 

2   !Jv       n 
_-iO   .   ,30  (H^)   -i(      °,,,oi'  ^-^^> 

d      U 
P 

The  expression  for  E[NCt)],  Eq.  IV-17,  can  now  be  substituted  into  the 
above  equation  to  obtain 


AP(0  _  ,.„  A  -   c  2  ^\ 


a  + 


L 


)  -^    [ _^.  ^  .^J  (IV-27) 


For  simplicity,  we  attribute  an  increase  in  the  pressure  drop  only 

to  the  increase  in  the  linear  velocity  of  flow;  this  is  equivalent  to 

2 
assuming  that  the  porosity  term,  [  (1-e) /e]    ,  in  Eq.  IV-27  remains  unchanged. 

Equation  IV-2  7 ,  therefore ,  reduces  to 


AP(t)  ^  ^_  AP^   f    g  +  B   _^  1  (IV-28) 


where  (-AP/L).  is  a  constant  and  reflects  the  Initial  pressure  drop 
per  unit  length  of  the  filter  bed. 

COMPARISON  OF  THE  MODEL  WITH  THE  EXPERIMENTAL  DATA  AND  DISCUSSION 

The  experimental  results  by  previous  researchers  (Eliassen,  1935; 
Ives,  1961;  Rimer,  1968;  Deb,  1969;  Huang,  1972)  have  been  analyzed  based 
on  the  present  model;  the  results  are  shown  in  Figures  IV-1  through  IV-11. 

The  values  of  a  and  g  in  the  model  equation,  Eq.  IV-28,  have  been 
determined  by  the  following  procedure.   The  preliminary  value  of  a  is 
first  determined  from  the  initial  slope  of  the  time  dependent  pressure 
drop  curve.  With  this  value  in  hand,  the  preliminary  value  of  0   is  estimated 
from  the  entire  set  of  the  data.   Finally,  the  values  of  a  and  6  are  recovered 
by  using  this  set  of  preliminary  values  as  a  starting  point  of  search  for 
the  minimum  of  sum  of  squares  residual  by  means  of  a  non-linear  optimizing 
technique  (see,  e.g.,  Marquardt,  1963).  The  initial  pressure  drop  per  unit 
bed  depth,  (-AP/Dj,,  which  is  also  contained  in  the  equation,  can  be  estimated 
by  the  Carman-Kozeny  equation  or,  more  accurately,  can  be  measured  at  the 
onset  of  the  experimental  run.   The  measured  values  of  the  initial  pressure 
drop  have  been  used  in  the  present  work. 

Deb  (1969)  conducted  his  experiment  by  allowing  turbid  water  to  flow 
through  a  deep  sand  bed  with  uniform  grain  size  of  0.647  mm  at  a  constant 

filtration  rate.   The  concentration  of  Fuller's  earth  particles  in  turbid 
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water  was  45  x  10  vol7o.   The  pressure  drop  were  measured  at  various  depths 

during  the  filtration  run  without  disturbing  the  normal  operation.   The 

present  model,  in  which  the  scouring  effect  is  negligible,  giving  rise  to 
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a  pure  birth  process,  appears  to  adequately  explain  the  increase  in  the 
pressure  drop  (see  Figure  IV-1) .   The  negligible  scouring  effect  seemed 
to  be  due  to  the  short  operating  time  as  well  as  the  diluteness  of  influent 
concentration,  for  which  the  deposit  accumulation  within  the  bed  was  small 
so  that  the  blockage  mechanism  became  far  more  dominant  than  the  scouring 
mechanism. 

A  rapid  filter  with  a  depth  of  60  cm  was  used  by  Eliassen  (1935)  to 
remove  suspended  hydrous  ferric  floe  particles  from  the  water.   The  pressure 
drop  reading  was  taken  every  ten  hours;  the  results  are  shown  in  Figure  IV-2. 
Each  experimental  run  lasted  sufficiently  long  to  retain  a  large  amount  of 
deposits  in  the  filter.  Under  such  conditions,  a  significant  scouring 
effect  probably  took  place - 

The  present  model  shows  a  good  fit  with  the  data  of  Huang  (1972), 
which  appear  in  Figures  IV-3  through  IV-5;  the  only  exception  is  run  B-5-III 
in  Figure  IV-5.  This  deviation  was  probably  caused  by  the  shallow  penetration 
of  suspended  particles  along  the  bed.   Huang  measured  the  concentration  of 
suspended  particles  at  different  depths  of  the  bed  under  various  flow  rates 
as  illustrated  in  Figure  IV-6.   Notice  that  when  the  flow  rate  was  4.9  m/h, 
the  filter  medium  beyond  a  depth  of  35  cm  did  not  contribute  to  the  further 
removal  of  suspended  particles,  rendering  the  lower  15.8  cm  of  the  bed 
unutilized.  Thus,  we  can  postulate  that  the  present  model  is  only  applicable 
to  the  first  35  cm  of  the  bed,  as  illustrated  in  Figure  IV-5. 

The  results  of  Rimer  (1968)  are  shown  in  Figure  IV-7  and  IV -8,  and 
those  of  Ives  (1961)  in  Figure  IV-9.  The  present  model  agrees  well  with 
the  data  of  Rimer  from  the  entire  bed  depth.   However,  the  model  can  describe 
Ives'  data  only  from  the  first  portion  of  the  bed.   This  is  indicative  of 
the  shallow  penetration  discussed  in  the  preceding  paragraph.   Ives  measured 
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in  detail  the  specific  deposit  as  a  function  of  time  and  bed  positions 
by  means  of  the  radioactive  algae  technique  and  determined  conclusively 
that  the  penetration  was  indeed  shallow. 

The  data  of  Deb  and  Rimer  show  that  the  operating  times  in  their 
experiments  were  roughly  equal  (see  Figures  IV-1  and  IV-7) .   However, 
the  pressure  drop  through  Deb's  filter  was  modeled  as  a  pure  birth 
process  but  that  through  Rimer's  filters  was  modeled  as  the  birth-death 
process,  for  which  the  scouring  effect  was  taken  into  account.  The 
need  to  incorporate  the  scouring  effect  in  the  model  for  Rimer's  data 
can  be  justified  in  the  light  of  measurement  of  effluent  concentration. 
In  Rimer's  filter,  the  effluent  concentration  continued  to  increase  with 
time  as  shown  in  Figure  IV-10.   This  is  indicative  of  the  scouring 
phenomenon;  recall  that  the  scouring  rate  increase  with  an  increase  in 
the  deposit  accumulation  as  can  be  seen  in  Eq.  IV-6.   In  contrast.  Deb's 
effluent  concentration  data  illustrated  in  Figure  IV-11  do  not  exhibit 
such  a  monotonically  increasing  trend. 

The  present  model  equation,  Eq.  IV-28,  contains  two  adjustable 
parameters,  a  and  3,  These  two  parameters  are  functions  of  many  factors, 
among  which  are  the  size  distributions  of  collector  grains  and  suspended 
particles,  properties  of  the  liquid  and  involved  surfaces,  filtration 
rate,  bed  porosity,  and  suspension  concentration. 

The  parameters,  a   and  6,  increase  with  an  increase  in  the  filtration 
rate  as  shown  in  Figure  IV-4.   A  plausible  reason  is  that  at  the  same 
influent  concentration,  the  filter  with  a  higher  flow  rate  is  likely  to 
retain  more  solid  particles  than  that  with  a  lower  flow  rate  within  the 
s  ame  t  ime  pe  r i  od . 
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The  parameters  are  functions  of  the  suspension  concentration  which 
decreases  along  the  bed  during  a  filtration  run.   Thus,  these  two  parameters 
might  also  decrease  with  the  bed  depth  as  indicated  in  Figure  IV-7.   In 
contrast,  a  decreases  monotonically  with  the  bed  depth  while  3  varies 
irregularly  as  indicated  in  Figure  IV-8,   Notice  that  the  bed  in  Figure 
IV-7  is  of  unisized  sand  while  the  bed  in  Figure  IV-8  is  of  multimedia 
(anthracite,  sand,  and  garnet).   It  appears  that  a  depends  strongly  on  the 
suspension  concentration  while  6  depends  weakly  on  the  suspension  concentra'^ 
tion;  the  latter  is  also  highly  dependent  on  the  type  of  bed  grains  and  bed 
structure,  i.e.,  the  spatial  distribution  of  these  grains. 

The  stochastic  birth-death  model  considered  here  is  assumed  to  be 
time-homogeneous,  in  which  both  parameters,  a  and  6,  are  not  functions  of 
time  because  of  the  constant  filtration  rate  and  suspension  concentration., 
i.e.,  the  constant  operating  conditions.   For  the  case  where  the  filtration 
rate  or  suspension  concentration  varies  with  time  during  a  filtration  run, 
the  birth-death  model  is  still  applicable,  but  both  parameters  will  be 
functions  of  time;  in  other  words,  the  process  described  by  the  model  is 
the  time-heterogeneous  process. 
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NOMENCLATURE 

d  =  diameter  of   Che   anthracite 

a 

d  =  average  diameter  of  the  grain  particles 
P 

d  =  diameter  of  the  sand 
s 

E[N(t)]  =  mean  value  of  the  random  variable  N(c) 

G(s,t)  =  probability  generating  function 

L  =  depth  of  the  bed 

ti  =  number  of  pores  clogged  per  unit  volume  of  the  bed 

n  =  total  number  of  open  pores  susceptible  to  blockage  per  unit 

volume  of  the  bed 

N(t)  =  a  random  variable  which  describes  the  number  of  blocked  pores 

per  unit  volume  of  the  bed  at  Che  moment  t 

p  (t)  =  probability  that  exactly  n  pores  are  blocked  at  the  moment  C 

(-AP/L)  =  pressure  drop  per  unic  length  of  the  filter 

2 

R  =  fraction  of  variation 

r  =  radius  of  pore 

s  =  variable  of  the  probability  generating  function 

t  =  time 

u  =  superficial  velocity 

V  .=  linear  velocity 

Var[N(t)]  =  variance  of  the  random  variable  N(t) 

z  =  height  along  the  filter  bed 

J  =  intensity  of  the  birch  cransicion 

V  =  intensity  of  the  death  Cransitlon 
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a      =  proportionality  constant  defined  In  Eq.  IV-5,  the  blockage 

constant 
Q  =  proportionality  constant  defined  in  Eq.  IV-6,  the  scouring 

constant 
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Figure  IV-1.   Fitting  the  present  model  to  Deb's  data  (1969): 

sand  size  (0.0647  cm),  bed  depth  (61  cm),  flow  rate  (4.9  m/h) , 

-4 
Fuller  s  earth  suspension  (45x10   vol%). 

Run  1;   position  0.2  (a=  0.261  h" 


position  0.4  (a  =  0.217  h 

position  0.6  (a  =  0.164  h 

position  0.8  (a  =  0.146  h" 

position  1.0  (a-  0.130  h~ 


S 


0.000   h 


■1 
,  -1 


R     =    0.992), 


0.000  h  ,  R  =  0.990), 
0.000  h"""",  R^  =  0.986)  , 
0.000  h"""",    R^  =    0.987). 
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Figure  IV-2.  Fitting  the  present  model  to  Eliassen's  data  (1935): 

sand  size  (0.061  cm),  bed  depth  (60  cm),  flow  rate  (4.9  m/h), 
hydrous  ferric  oxide  floe  suspension  (concentration ,  7.3x10 


wt%). 
Run  6; 


a   =  0.0252  h 


R  =  0.996. 
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Figure  IV-3.  Fitting  the  present  model  to  Huang's  data  (1972): 

bed  depth  (60.96  cm),  flow  rate  (9.8  m/h),  waste  water  sus- 
pension (12.5  mg/£). 


Runs;  C-3-I   (a  =  0.645  h 
C-3-II   (a  =  0.578  h' 


=  0.0146  h  "'",  R  -  0.997), 


C-3-III  (a  =  0.382  h 
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Figure  IV-4.   Fitting  the  present  model  to  Huang's  data  (1972): 

bed  depth  (60.96  cm:  anthracite,  30.48  cm,'  d  =  0.184  cm;  sand, 
30.48  cm,  d  =  0.055  cm),  waste  water  suspension  (12.5  mg/i) . 


Runs;   F-l-I   (a  =  0.122  h 
F-l-II   (a  =  0.100  h" 


F-l-III  (a  =  0.089  h 


=  0.0054  h  ^,  R^  =  0.999), 


=  0.0045  h"""",  R^  =  0.992), 


=  0.0040  h  •^,  R^  =  0.993). 
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Figure  IV-5.  Fitting  the  present  model  to  Huang's  data  (1972): 

sand  size  (0.092  cm),  bed  depth  (50.8  era),  waste  water  suspension 
(12.5  mg/£). 


Runs;  B-5-I    (a 
B-5-II   (a 


0.787  h 


-5-III  (a  =  0.698  h 


R  =  0.996), 
R^  =  0.997), 
R^  =  0.999). 
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Figure  IV-6.   Extent  of  suspended  particles  removal  per  cm  depth  of 
a  unisized  sand  filter: 
sand  size,  0,092  cm;  filtration  time,  1  h; 


Runs;  B-5-I,  IX,  III  (Huang,  1972). 
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Figure  IV-7.   Fitting  the  present  model  to  Rimer's  data  (1968): 

sand  size  (0.046  cm),  bed  depth  (60.96  cm),  flow  race  (7.3  m/h) , 
Feci,  suspension  (5  mg/H,  pH  =  8.3). 

-1  -1    2 

Run  VII-2;  position   4.24  cm  (a  =  1.035  h   ,  B  =  0.069  h   ,  R  =  0.986), 

-1  -1    2 

position  13.47  cm  (a  =  0.718  h   ,  3  =  0.066  h   ,  R  =  0.981), 

-1  -1    2 

position  23.70  cm  (a  =  0.511  h   ,  B  =  0.063  h   ,  R  =  0.997), 

-1  -1    2 

position  60.96  cm  (a  =  0.347  h  ,  B  =  0.060  h   ,  R  =  0.998). 
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Figure  IV-8.  Fitting  the  present  model  to  Rimer's  data  (1968): 

bed  depth  (60.96  cm:  anthracite,  20.32  cm,  d  =  0.100  cm;  aand, 
d  =  0.071  cm;  garnet,  20.32  cm,  d  =  0.059  cm),  flow  rate  (7.3 


position  34.72  cm  (a  =  0.319  h 


20.32  cm, 

m/h). 
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Figure  IV-9.  Fitting  the  present  model  to  Ives'  data  (1961): 

sand  size  (0.05AA  cm),  bed  depth  (20  cm),  flow  rate  (A. 9  m/h) , 
chlorella  algae  suspension  (135x10   vol%). 

Run  IB;  position  9.2  cm  (a  =  0.359  h"  ,  3  =  0.102  h~  ,  r'  =  0.986). 
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Figure  IV-10.  Ratio  of  the  effluent  to  influent  ferric  iron  concnetrations, 
C  /Cq,  as  a  function  of  time:  Run  VII-2  (Rimer,  1968). 
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Figure  IV-11.   Isoconcentration  ratio  lines:  Run  1  (Deb,  1969). 


APPENDIX  A:   AN  ALTERNATIVE  DERIVATION  OF  THE  MASTER  EQUATIONS 

The  Master  equations,  Eqs.  IV-3  and  IV-4,  in  the  text  can  also  be  derived 
through  the  consideration  of  a  time  homogeneous  Markov  process.   Let 
p..(T,t)  denote  the  transition  probability  that  the  system  in  the  state 
of  i  pores  blocked  at  time  t  will  change  to  the  state  of  j  pores  blocked 
at  time  t;  i,J  =  0,1,2, 
time  invertal  Ct,t+At) 


p..(t,t+4t)  =  v.. At  +  o(At)  ,  ij'J;  i,j  =  0,1,2,.. .,n      (A-1) 


1  -  p..(t,t+At)  =  -v..At  +  o(At)  ,  i  =  0,1,2 n         (A-2) 

where  v. .  and  v. .  are  the  intensities  of  transition, 
ij      li- 
lt can  be  shown  that  the  transition  probabilities,  {p..(T,t)} 

satisfy  the  following  Kolmogorov  forward  differential  equations  (Chiang, 

1980), 


ft   "ij^^'^^  =  ^}^^^''"''   'M  •  ''^    '   0.1. 2. ■■-."0'        (^-3) 
with  the  initial  condition 

fl,       i  =  J 
p.j(T,T)  =  (A-4) 

[o,       i  )*  j 

Let  p..(t)  denote  the  transition  probability  from  time  0  to  time  t. 
In  matrix  notation,  Eqs.  A-3  and  A-4  can  be  written,  respectively,  as 


(A-5) 
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P(0)  =   I  '   Identity  matrix 


(A-6) 


P(t) 


P  „(t:)   P   ,(t) 


'<-' 
W 


P  „  (t) 


(A-7) 


"00 
^0 


lV 


01 

^1 


"On, 


(A-8) 


In  the  present  model,  it  is  assumed  that  given  NCt)  =  n 
1.   The  transition  probability  that  more  than  one  pore  will  be  blocked 
or  scoured  is  o(At) ,  i.e.. 


=  0 


|j  -  i|  >  2;     i,j  =  0,1, 


n..  (A-9) 

ij      '     1^     I  -   ■      ■■'  0 

2.   The  transition  probability  that  an  open  pore  will  be  blocked  during 

the  interval  (t,  t+At)  is  X  At  +  0(At),  i.e., 
n 

\      ,  i  =  n  ;      n  =  0, 1 , 2 , . . . , n.-l  (A-10) 

n  0 

j  =  n  +  1 


IJ 
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3.   the  transition  probability  that  a  blocked  pore  will  be  scoured 
during  the  interval  (t,  t+At)  is  u  At  +  o(At) ,  i.e., 

v,  =  u,     i  =  n;      n=  1,2 n.;  (A-11) 

IJ     n  '  0 

j  =  n  -  1 
A.   The  transition  probability  of  the  no  change  in  the  interval  (t,  t+At) 
is  II  -  X  At-  p„At  -  o(At)l,  i.e.. 


-(A   +  M  )  ,     i  =  n;      n  =  0,1,2, 


(A-12) 


Then  Eq.  A-8  reduces  Co 


^0 

0 

0 

-(A^+U,) 

h 

0 

"2 

-(^2+^2) 

A 

)        0  0  0        .    .    .        -V 

n 

Let  p  (t)  be  the  probability  that  exactly  n  pores  are  blocked 
at  time  t,  and  fCt)  the  vector  of  probabilties ,  i.e., 


(A-13) 


Jl(t)  =  Ip.(t)  p  (c)  p  (t)  ...  p   (t)J 
u      J.      i  n^ 


(A-14) 


It  can  be  seen  that 

^(t)  =^(0)P(t)  (A-15) 

Multiplying  both  sides  of  Eq.  A-5  with  P  (0)  ,  we  obtain 


(A-16) 
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-^  P  (c)  =  lPCt)U  C16a) 

dt  — 

Expanding  Eq.  A-16a,  we  obtain 

dp  (t) 


dPott) 


-A  P„(t)  +  v.? At)  (A-18) 


dt    "0  0''   "^1  r 

which  are  the  same  as  Eqs.  IV-3  and  IV-4,  in  the  text. 
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APPENDIX  B:   SOLUTION  OF  THE  PROBABILITV  GENERATING  FUNCTION  OF  THE 
MASTER  EQUATIONS 

The  system  of  differential  equations,  Eqs.  IV-7  and  IV-8,  in  the  text, 

describe  the  linear  stochastic  birth-death  process  as 


dp  (t) 

—j^     =  c.[n   -  (n  -  l)]p_^_^(t)  -  [a(n^  -  n)  +  Bn]p^(t) 

n  >  1  (B-1) 


dp^(t) 

The  mean  number  of  the  blocked  pores  at  a  given  moment  c  is  to  be 
determined;  it  is 

n 
E[N(t))  =   E°  np  (t)  (B-3) 

n=0 

To  evaluate  this  value,  we  use  the  method  of  probability  generating 
function  which  is  defined  as 

G(s,t)  =   "°  s"p  (t)  (B-4) 

n=0    " 

Multiplying  both  sides  of  Eq .  B-1  by  the  respective  s  's  and  Eq .  B-2 
by  s  ,  the  following  expression  is  obtained; 


"U  n  "Pn^^^    ,„  ,  ,     „,       2„  ^O  „  n-1 


z"  s"  ^V-  =  [6  +  (d  -  B)s  -  as"^][  l"  ns"""p  (t)] 
n=0      ^^  n=0 


+  [an^(s  -  1)1[  1:°  s"p  (t)l  (B-5) 
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Ms^  =  J  s"— 2— -  (B-6) 


aG(s,t)    ^0   n-l   ,  ,  ,„  7, 

— „  '  '    =  Z  ns   p  (t),  \o-i) 

3s        „       n 
n=0 


Equation  B-5  reduces  to 


aG( 

3t 


^i^^  =  -^^If^  [e  +  Cc-  B)s  -  cxs^]  +  GCs,t)[c<ng(s  -1)1   (B- 
The  Initial  and  boundary  conditions  can  be  expressed  as 


G(s,0)  -  1  (B-9) 

G(l,t)  =  1  (B-10) 

Assume  that  the  solution  of  Eq.  B-8  take  the  form  of 

n 
G(s,t)  =  [f(3)  +  g(s,t)]  °  (B-11) 

Substituting  this  into  Eq.  B-8  yeilds 

M|^=I,.(„-3)s-as^lI^.^^^] 

+  [a(s-  l)][f(s)  +g(s,t)l  (B-12) 

Let  f(s)  and  g<s,t)  satisfy 

[B  +  (a  -  S)s  -  as^]-^  +  [o(s  -  l)]f(s)  =  0  (B-13) 

ds 


'^^^   =  [B  +  (a  -  B)s  -  as^]-^^^-^  +  l°(=  -l)]g(s,t)     (B-14) 
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From  Eqs.  B-9  and  B-11,  we  obtain 

g(s,0)  =  1  -  f(s)  (B-15) 

From  Eqs.  B-10  and  B-11,  we  have 

f(l)   =  1  (B-16) 

gCl,t)  =  0  (B-17) 

Subject  to  these  boundary  and  initial  conditions,  Eqs.  B-13  and  B-14 

are  solved  ,  respectively,  as 


£(s)  =  ^ 


(B-18) 


,     ,,  -{a   +  B)t 
g(s,t)  =  ^^,r^ (B  19) 

Therefore,  the  solution  of  Eq.  B-18  becomes 

G(s,t)  =  ;(«s+  B)  -  a(s  -  De''"-^  "  Y°  (B-20) 

which  is  the  same  as  Eq.  I¥-15  in  the  text. 


CHAPTER  V   BIRTH-DEATH  MODELING  OF  DEEP  BED  FILTRATION: 
Sectional  Analysis 

In  the  stochastic  modeling  presented  in  Chapter  IV,  or  more 
specifically  birth-death  modeling,  it  has  been  assumed  that  the  intensity 
of  birth  or  the  blockage  intensity,  A  ,  and  the  intensity  of  death  or  the 
scouring  Intensity,  p  ,  are  linearly  proportional  to  the  number  of  open 
pores  and  the  number  of  blocked  pores  in  the  bed,  respectively  (Litwiniszyn, 
1968a;  1958b;  1969;  Hsu  and  Fan,  1984).   This  assumption  of  linearity  gives 
rise  to 

^n  '   °'"0  "  "'-         n  =  0,  1,  2.  ....  ng  (v-1) 

U   =  6n;  n  =  0,  1,  2,  . . . ,  n^  (V-2) 

where  n  is  the  total  number  of  open  pores  susceptible  to  blockage  in  a 
unit  volume  of  the  bed,  a  is  the  blockage  constant,  and  3  is  the  scouring 
constant.   Equation  V-1  is  based  on  the  postulation  that  as  the  number  of 
open  pores  susceptible  to  blockage  becomes  larger,  the  probability  of  an 
open  pores  being  blocked  will  be  proportionally  greater,  and  Eq .  V-2  on 
the  postulation  that  as  the  number  of  blocked  pores  becomes  larger,  the 
probability  of  a  blocked  pore  being  scoured  will  be  proportionally  greater. 
The  applicability  of  this  modeling  to  the  available  experimental  data  has 
been  amply  demonstrated  in  Chapter  IV. 

It  may  be  noted,  however,  that  the  linear  velocity  of  flow  through  the 
bed,  V,  will  increase  when  the  number  of  open  pores,  (n„-n) ,  decreases;  in 
fact,  V  is  inversely  proportional  to  (n^-n) .   Furthermore,  solid  particles, 
suspended  in  faster  moving  fluid,  will  probably  be  less  likely  to  be  trapped 
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in  the  bed.   Thus,  we  may  also  postulate  that 

'\  °  V  "  ^"O  -  "^  (V-3) 

Because  of  two  effects  of  the  number  of  open  pores  on  A  ,  as  characterized 
by  Eqs.  V-1  and  V-3,  it  is  highly  probable  that  X   is  proportional  to  the 
square  of  the  number  of  open  pores  under  some  conditions.   In  other  words, 
under  this  second-order  approximation,  A   assumes  the  form 

Similarly,  y  assumes  the  form 

%  '   8"^  (V-5) 

MODEL 

A   second-order,    birth-death  model   represented   by  Eqs.    V-4  and  V-5   is 
employed    to  describe  both  blockage   of   pores   by   suspended   particles   and 
scouring  of  deposited   particles.      This  model,    coupled  with   the  Carman- 
Kozeny  equation,    simulate    the   pressure  drop  dynamics  during  a  deep  bed 
filtration. 

Master  Equation  and  One-step  Operator 

The   number  of  blocked    pores    in  a  unit   volume   of   the  bed   at    time    t, 
N(t) ,    is   taken  as    the   random  variable    (Litwiniszyn,    1966)  ;    a   specific 
value   of   N(t)   will  be  denoted    by  n.      For   the   birth-death   process,    repre- 
senting  the  blockage-scouring   process    of   pores    in   the    filter  bed,    the 
following   is    assumed    (Chiang,    1980). 

1.      The   conditional    (or   transition)    probability   that   a   birth   event, 
i.e.,    the  blockage   of  an  open   pore,   will    occur  during    the    time    interval 
(t,    t+At)    is    X   At  +  o(At),   where    A      is   a    function  of  n,    (e.g.,    given  by 
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Eq.  V-1  or  V-4) . 

2.  The  conditional  probability  that  a  death  event,  i.e.,  the  scouring 
of  a  blocked  pore,  will  occur  during  the  time  interval  (t,  t+At)  is  y^At  + 
o(At),  where  y  is  also  a  function  of  n,  (e.g.,  given  by  Eq.  V-2  or  V-5). 

3.  The  conditional  probability  that  more  than  one  event  will  occur 
in  this  time  interval  is  o(At),  which  signifies  that 

.   o(At)  _ 

4.  The  conditional  probability  that  no  change  occurs  during  the  time 

interval  (t,  t+At)  is 

[1  -  A  At  -  p  At  +  o(At)] 

The  relationship  among  these  four  conditional  probabilities  can  be  diagrammed 
as  shown  in  Figure  V-1.   By  denoting 

p  (t)  =  Pr{N(c)  =  n};    n  =  0,  1,  2 n^ 

as  the  probability  that  n  pores  are  blocked  at  time  t,  the  probability 
balance  around  N(t)  =  n  yields  (Chiang,  1980) 


^  =  K-1   Vl(^>  ^  %+l  "n+l^^^  -  ^^n  ^  ^nK^^'-' 


n  =  1,  2,  (n.  -  1),  (V-6) 


^^  =  -XQPo(t)  +  u^P^(t)  (V-7) 


dp   (t) 

-^  -   '^n^-l  ^„-l'^'  -  \  \^^'  (V-8) 
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Letting 


Eqs.  V-6  through  V-8  can  be  combined  as 
dp  (t) 


dt 

n  =  0,  1,  2 n  (V-9) 

0 

Equation  V-9  is  the  so-called  probability  balance  equation  or  Master 
equation  (see,  e.g.,  Van  Kampen,  1981). 

The  one-step  operator, E ,  is  defined  by  its  effect  on  an  arbitrary 
function  f(n)  as  (see,  e.g..  Van  Kampen,  1981) 

(V-10) 

With  the  aid  of  this  operator,  the  Master  equation  for  the  birth-death 
process,  Eq.  V-9,  reduces  to 

^^^^     '      (E  -  l)y^p(n,t)  +  (]e'^  -  l)A^p(n,t)        ^^.^i) 

where  p  (t)  has  been  rewritten  as  p(n,t)  for  convenience. 

Expansion  of  Master  Equation 

A   and  u  ,  defined  by  Eqs.  V-4  and  V-5,  respectively,  are  nonlinear, 
n      n 

and  thus,  the  analytical  solution  of  Eq.  V-11  is  extremely  difficult,  if 
not  impossible.   It  has  been  shown,  however,  that  this  difficulty  can  be 
circumvented  by  resorting  to  approximation  based  on  the  Taylor's  expansion 
(see,  e.g.,  Van  Kampen,  1981). 


For   linear   X     and    n     defined   by  Eqs.   V-1  and  V-2,    respectively, 

n      n 

the  mean  of  N(t),  <n>,  obtainable  from  the  Master  equation,  Eq.  V-11, 
is  also  obtainable  from  the  corresponding  deterministic  governing 
equation  given  below  (also  see  APPENDIX  A) . 

^^i^  =  C.[njj-S(t))  -  Bn(t)  (V-12) 

where  n(t)  represents  the  deterministic  value  of  n  at  the  moment  t. 
Accordingly,  the  corresponding  deterministic  governing  equation  with 
the  second-order  approximation  becomes 

dn(t)  _   ,     -,  ,,2 


^^     aln^  -  n(t))   -  Sn'Cc)  (V-13) 

Letting 

Eq.  V-13  reduces  to 

4^=  an„[l-Kt)]^  -  BnyCt)  (V-14) 

at  U  u 

Knowing  that  n  has  a  maximum  value  of  n^,  we  expect  that  n  consists 

1/2 
of  the  macroscopic  tenn,  n  ^i  (t)  ,  and  the  fluctuating  term  of  order  n-   , 

that  is,  p(n,t)  will  have  a  sharp  peak  located  roughly  at  n^^iji  (t)  and  a 

1/2 
width  of  order  n-    (see  Figure  V-2).   Hence,  we  set 

1/2 
n  =  nQ$(t)  +  Hq   C  .  (V-15) 

where  ^  is  a  new  variable.  Notice  that  i|j(t)  satisfies  the  deterministic 
governing  equation  as  will  be  elaborated  later.   Accordingly,  the  proba- 
bility distribution  p  of  n  now  becomes  the  probability  distribution  tt 
of  i;.   Thus, 
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p(n,c)fln  =  Tl(C,t)A; 


This  in  conjunction  with  Eq.  V-15  yields 

iT(C.t)  =  ny%(n,t)  (V-16) 


The  time  derivative  in  Eq.  V-11  is  taken  with  constant  n;  This  means 
that  this  derivative  is  taken  in  the  5-t  plane  along  the  direction  given 
by  (see  Eq.  V-15) 

d£  ^  _  1/2  d$ 
dc     "o   dt 


1/2  3p(n,t)  _  OiT    1/2  d^  3£  ,,  ,,, 

"o     ~T^       '^  37  "  "o     dt  a;  ^      " 

-1/2 
The  operator  E  changes  n  into  n+1 ,  and  therefore,  ^  into  j+Hq   ,  so  that 

(see  APPENDIX  B) 

-1/2  3    1   -1  3^  _^ 
^  =  ^  +  "0    3^^2"0  ;T+  •••  (V-18) 

and 

,r.-l    ,    -1/2  3^1  -1  3^  _^ 

^        =  1  -  "O    i?""  2  "O  ~2  ""  •••  (V-19) 

3^ 

Substituting  Eqs.  V-16  through  V-19  into  the  Master  equation,  Eq.  V-11, 

-1/2 
tactor  n^ 

3-n  _   1/2  dj  311^ 
3t  '   "O   dt  3; 


dc. 


.an^(-n-l/2  A  +  in-l^-  ...)(!- 0  -  n'^'^^)  ".  (v-20) 
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By  absorbing  the  factor  n„  into  the  time  variable  as 


Eq.  V-14  is  simplified  as 


(V-21) 


-j^  =  a(l-(j>)  -  B  $  (V.22) 


The  transformation  given  by  Eq.  V-21  together  with  the  truncation  of  the 
terms  after  the  second  derivative  for  large  n-  gives 


3tt    1/2  d(t  2IL 
3t  "  "o   dx  3? 

„   ,  -1/2  3^1-13^,,*,   „-l/2  2 
=  2"0^"0    ^^  2"0  -^^    ^*.+  "O   ^^   " 

Expanding   the   right-hand   side   of  Eq.    V-23   and   collecting   the   resultant 
terms    of  order   of  n_        and   n_    separately,   we  have 

^=  a(l-0)^  -   6$^  (V-24) 


1^  =  2[3$  -aCl-*)]  -^  C^TT)   +1  13*^  +  cLd-'P)^]  -^  TT  (v-25) 

That  Eq.  V-24  is  the  same  as  the  deterministic  governing  equation,  Eq.  V-22, 
satisfies  our  postulation  that  4'  is  governed  by  the  deterministic  governing 

equation.   Equation  V-25  suffices  for  determining  the  first  and  second 

2 
moments  of  c; ,  <^>  and  <<;  >  ,  respectively,  which  contain  the  significant 

information  on  the  random  variable.   Upon  transformation  of  Eq.  V-25,  we 
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obtain  (see  APPENDIX  C) 

■^=  [2a(l-4>)  -  2B*]  <C>  (V-26) 


Since 


1/2 


as  indicated  by  Eq,  V-15,  we  have 


<n>  =  n  ij)  +  n^''^<?>  (V-28) 


(V-29) 

2 
where  <n>  is  the  first  moment  of  n,  <n  >  is  the  second  moment  of  n,  and 

2 
<<n  >>  is  the  variance  of  n. 

For  a  pure  birth  process  where  g  =  0,  Eq.  V-22  reduces  to 

d*  2 

-^  =  a(l-*)  (V-30) 

Solving  this  equation,  subject  f  the  initial  condition 

4)  =  0   at   T  =  0,  (V-31) 

gives 
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Substituting  this  expression  into  Eq.  V-26  and  setting  g  =  0  in  the 
resultant  expression,  we  obtain 

m^  =  _2£L_  <,>  (V-33) 

3t    1  +  ax  ^ 


The  general  solution  of  this  equation  is 


(V-34) 


To  satisfy  the  initial  condition 

<;>  =  0   at   T  =  0,  (V-35) 

the  constant  c  must  be  zero;  hence  (see  APPENDIX  D) , 

<;;>  =  0  (V-36) 

Substitution  of  Eqs .   V-32  and  V-35   into  Eq.    V-28   gives 


an  T 
<">  =  T-rT7  (V-37) 


Knowing  that 


as  indicated  by  Eq.  V-21  and  letting 

an^  =  k  (V-38) 

Eq.  V-37  reduces  to 

"">   =  "n  T-^  (V-39) 
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The  expression  is  a  reasonable  approximation  to  <n>  for  a  large  n^,  where 
the  terms  after  the  second  derivative  in  Eq.  V-20  are  negligible,  and 
therefore,  it  is  valid  for  the  deep  bed  filtration  process.   The  variance 
of  n  can  also  be  obtained  by  solving  Eqs .  V-27  and  V-29.   For  a  sufficiently 
large  n_,  the  variar 
filtration  process. 

Constant  Flow  Filtration 

Denoting  the  linear  velocity  of  the  suspension  flowing  through  the 
pores  at  time  t  as  v(t)  and  the  linear  velocity  at  the  onset  of  filtration 
we   have 

under  the  assumption  that  all  pores  have  the  same  radius. 

For  the  second-order,  pure  birth  process,  substitution  of  Eq.  V-39 
into  Eq.  V-40  yields 

(V-41) 

Under  the  simplifying  assumption  that  an  increase  in  the  pressure  drop 
can  be  attributed  only  to  the  increase  in  the  linear  velocity  of  flow, 
substituting  Eq.  V-41  into  the  Cannan-Kozeny  equation, 


AP(t)  _  ,^„  ,l-€,^   tiv(t) 


,    =  150  i.-^^-^)      ii^^      ,  (V-42) 

^  ^     d 

P 


yields 
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2  uv 
^  ^      d 


(_AP)   (l  +  kt)  (V-43) 


0 


where  [-(iP/L).]  represents  Che  initial  pressure  drop  per  unit  length  of 
the  bed.   For  the  linear,  birth-death  modeling  presented  in  Chapter  IV, 
the  pressure  drop  through  the  filter  bed  has  been  shown  to  change  as 
(also  see  APPE^!DIX  E) 


AF(t)  _  ,  AP,  j    a  +  e,  I 


(V-44) 


ANALYSIS  OF  AVAILABLE  EXPERIMENTAL  DATA  AND  DISCUSSION 

Many  factors  may  affect  the  blockage  and  scouring  constants,  a  and  fi; 
they  include  the  properties  and  concentration  of  suspension,  filtration 
rate  and  bed  characteristics.  Within  certain  ranges,  it  is  possible  to 
manipulate  independently  all  these  factors  except  the  suspension  concen- 
tration which  will  decrease  along  the  bed  due  to  the  trapping  of  solid 
particles.   It  is  entirely  conceivable,  therefore,  that  under  some  situations, 
a  lumped  model  treating  the  entire  bed  as  a  single  compartment  will  fail  to 
represent  the  deep  bed  filtration.   Thus,  it  may  be  desirable  to  divide 
the  bed  along  its  axis  into  more  than  one  compartment  and  to  measure  the 
increase  in  the  pressure  drop  through  each  compartment  under  such  situ- 
ation.  The  experimental  results  of  previous  researchers  (Eliassen,  1935; 
Ives,  1961;  Rimer,  1968;  Deb,  1969;  Huang,  1972)  have  been  analyzed  in  this 
light  based  on  the  present  models;  the  results  are  shown  in  Figures  V-3 
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through  V-9. 

In  Eliassen's  experiments,  a  rapid  filter  with  a  depth  of  60  cm' 
was  used  to  remove  suspended  particles  of  hydrous  ferric  oxide  floe  from 
suspension.  The  pressure  readings  were  taken  every  ten  hours  from  different 
bed  heights  (4.27  cm,  11.89  cm,  42.36  cm,  and  60  cm).   For  convenience, 
therefore,  the  bed  has  been  divided  into  four  compartments,  and  the  pressure 
drop  through  each  compartment  and  its  changes  have  been  simulated  according 
to  the  present  models..   The  results  of  runs  6  and  8  are  shown  in  Figures 
V-3  and  V-4,  respectively.  These  figures  show  that  the  change  in  the 
pressure  drop  through  the  first  compartment  (0  ~~   4,27  cm)  obeys  essentially 
the  second-order,  pure  birth  process  (Eq.  V-43) ,  that  through  the  second 
compartment  (4.27  —  11.89  cm)  reflects  the  linear,  birth-death  behavior 
(Eq.  V-44) ,  and  that  through  the  third  compartment  (11,89  -   42.36  cm)  follows 
the  linear,  pure  birth  model  (Eq.  V-4.4  with  0=0).  These  figures  also  show 
that  the  change  in  the  total  pressure  drop  through  these  three  compartments 
(0 -*  42.36  cm)  can  be  satisfactorily  simulated  by  the  lumped,  linear,  birth- 
death  model  for  which  X      and//   are  given  by  Eqs .  V-1  and  V-2,  respectively. 
Notice  that  the  pressure  drop  through  the  last  compartment  (42.36'^  60  cm) 
remained  essentially  constant  during  each  filtration  run.  This  implies 
that  the  last  compartment  was  not  utilized  because  of  the  relatively  shallow 
penetration  of  the  suspended  particles;  naturally,  the  present  models  are 
not  applicable  to  the  last  compartment  where  no  filtration  occured. 

The  data  of  Huang  (1972)  are  shown  in  Figure  V-5  and  those  of  Rimer 
are  shown  in  Figures  V-6  and  V-7;  they  both  measured  the  pressure  drops  at 
different  bed  heights.   These  figures  demonstrate  that  compartmentalization 
of  each  filter  leads  to  successful  simulation  of  its  dynamic  behavior  by 
the  present  models.   The  change  in  the  pressure  drop  through  the  compartments 
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around  the  entrance  section  follows  the  linear,  birth-death  model  and 
that  through  the  compartments  away  from  the  entrance  section  reflects 
the  linear,  pure  birth  characteristics.  These  figures  also  show  that 
the  change  in  the  overall  pressure  drop  exhibits  the  linear,  birth-death 
behavior.   Notice  that  in  Figure  V-6  the  pressure  drop  remained  constant 
through  the  last  compartment  of  Rimer's  filter  (44.45  ~  69.96  cm),  thus 
indicating  that  this  portion  was  not  utilized. 

The  results  obtained  by  Deb  (1969)  and  Ives  (1961)  are  plotted  in 
Figure  V-8  and  V-9,  respectively.   The  penetration  through  their  filters 
was  indeed  shallow.   It  can  be  seen  in  these  figures  that  in  each  case 
the  pressure  drop  remained  constant  through  more  than  507,  of  the  bed  away 
from  the  entrance  section.  However,  in  Deb's  filter,  the  pressure  drop 
through  the  remainder  of  the  bed  exhibited  the  linear,  pure  birth  behavior, 
and  in  Ives'  filter,  the  pressure  drop  through  the  remainder  of  the  filter 
reflected  the  linear,  birth-death  characteristics. 

In  this  work  both  sectional  analysis  and  lumped  analysis  have  been 
performed  and  the  results  are  compared.  The  sectional  analysis  gives  an 
insight  into  the  performance  of  the  deep  bed  filter  but  the  lumped  analysis 
provides  a  convenient  way  to  simulate  dynamics  of  the  overall  pressure  drop 
variation. 

According  to  Ellassen's  data  plotted  in  Figures  V-3  and  V-4,  the 
second-order,  pure  birth  behavior  is  found  to  occur  in  the  compartment 
around  the  entrance  section  (0  --  4.27  cm)  of  his  filter  bed.  Although  only 
77o  of  the  bed  exhibits  this  behavior,  its  contribution  to  the  increase  in 
the  overall  pressure  drop  is  almost  307o 

The  blockage  and  scouring  constants,  a  and  fi,   are  functions  of  the 
suspension  concentration  which  decreases  along  the  bed  during  a  filtration 
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run.   Thus,  If  the  bed  characteristics  are  uniform  throughout  the  bed, 
these  two  constants  are  expected  to  decrease  with  the  bed  depth  as  indicated 
in  Figures  V-5  and  V-6.  While  O  decreases  momotonically  with  bed  depth,  0 
varies  irregularly  as  seen  In  Figure  V-7.   Notice  that  the  bed  in  Figure 
V-7  is  of  multimedia  (anthracite,  sand,  and  garnet).   It  appears  that  a 
depends  strongly  on  the  suspension  concentration  while  fi   depends  weakly 
on  it;  however,  the  latter  is  highly  dependent  on  the  type  of  bed  grains 
and  the  bed  structure,  i.e.,  the  spatial  distribution  of  these  grains. 

As  mentioned  earlier,  the  pressure  drop  dynamics  in  the  section  of 
the  filter  beds  away  from  the  entrance  section  tend  to  obey  the  linear,  pure 
birth  process  as  shown  in  Figures  V-3  through  V-7.  A  plausible  interpretation 
is  that  in  these  sections,  the  suspension  concentrations  beenee  sufficiently 
dillute  so  that  the  blockage  mechanism  predominates  over  the  scouring 

mechanism.   This  is  equivalent  to  the  situation  seen  in  Figure  V-8;  the 

-4 
influent  suspension  concentration  was  only  45  x  10   volume7o  so  that  the 

scouring  effect  was  negligible,  which  gave  rise  to  the  linear,  pure  birth 

process . 
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NOMENCLATURE 

d  =  diameter  of  Che  anthracite 

a 

d  =  average  diameter  of  the  grain  particles 
P 

d  ■   diameter  of  the  sand 

s 

L  =   depth  of  the  bed 

n  =   number  of  pores  clogged  per  unit  volume  of  the  bed 


"0 


=  total  number  of  open  pores  susceptible  to  blockage  per  unit 
volume  of  the  bed 


N(t)  =  a  random  variable  describing  the  number  of  blocked 
pores  per  unit  volume  of  the  bed  at  the  moment  t 

< n>  =  mean  or  f i rs t  momen t  o f  n 

2 

<n  >  =  second  moment  of  n 

2  .       . 

<<n  >>  =  variance  of  n 

p  (t)  =  probability  that  exactly  n  pores  are  blocked  at  the  moment  t 

n 

pCn.t)  =  same  as  p  (t) 

(_  4£.)  =  pressure  drop  per  unit  length  of  the  filter 

r  =  radius  of  pore 

t  =  time 

u  =  superficial  velocity 

V  =  linear  velocity 

E  =  symbol  of  one-step  operator 

r  =  variable  representing  the  fluctuating  terra 

Ti(^,t)  =  probability  density  function 

X  =  intensity  of  the  birth  transition 
n 

n  =  intensity  of  the  death  transition 
n 

a  =  proportionality  constant  defined  in  Eq .  V-1  or  V-4,  the 
blockage  constant 

0  =  proportionality  constant  defined  in  Eq .  V-2  or  V-5,  the 
scouring  constant 

T  =  newvariabledefinedinEq.  V-21 
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Figure  V-1.   Birth-death  process  and  its  conditional  probabilities. 


V-17 


P(n,o) 


P(n,t) 


P"(n,t) 


n.0(o) 


Ho  </•(!) 


n.  V    (t) 


Figure  V-2.   Evolution  of  the  probability  density. 
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Figure  V-3.   Fitting  the  present  model  to  Eliassen's  data  (1935): 

sand  size  (0.061  cm),  bed  depth  (60  cm),  flow  rate  (^,9  m/h) ,  hydrous 
ferric  oxide  floe  suspension  (concentration,  7.3x10"  wt%)  . 

Run  6;  O  :  first  compartment  (0  '^   A.  27  cm),  k=0.123  h   , 

A  :  second  compartment  (4.2  7  i^  11.89  cm),  a=0.0A7  h"  ,  B=0.004  h~ 


O  :  third  compartment  (11.89  '^   42.36  cm),  a=0.011  h 


B=0, 


D 


lumped  model  (0  "V  42.36  cm),  a=0.032  h 
fourth  compartment  (42.36  ^   60  cm). 


6=0.0045  h 
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Figure  V-4.   Fitting  the  present  model  to  Eliassen's  data  (1935): 

sand  size  (0,061  cm),  bed  depth  (60  cm),  flow  rate  (4^9  m/h) ,  hydrous 
ferric  oxide  floe  suspension  (concentration,  7.3x10   wt%) . 

Run  8;  O  :  first  compartment  (0  "^   4.27  cm),  k=0.14  h   , 

A  :  second  compartment  (4.27  '^   11.89  cm),  a»^O.D373  h"  ,  8=0.0043  h""^ , 

O:  third  compartment  (11.89  "^   42.36  cm),  a=0.01  h"  ,  6=0, 

D  :  lumped  model  (0  "^   42.36  cm),  a=0.032  h"  ,  6=0.0058  h"  , 

0  :  fourth  compartment  (42.36  "^   60  cm). 
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Figure  V-5.   Fitting  the  present  model  to  Huang's  data  (1972); 

sand  size  (0.065  cm),  bed  depth  (33.35  cm),  flow  rate(4.9  m/h) ,  waste 

water  suspension  (12.5  tng/z  )  . 

first  compartment  (0  '^  12.7  cm),  a=0.429  h"  ,  6=0.02  h'  , 
second  compartment  (12.7  "^   33.35  cm),  a=0.232  h   ,  8=^0, 


Run  A-2-II;  O 
A 
O 


lumped  model  (0  ■^  33.35  cm),  a=0.368  h 


,  -1 


B=0.016  h 
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Figure  V-6.   Fitting  the  present  model  to  Rimer's  data  (1968): 

sand  size  (0.046  cm),  bed  depth  (60.96  cm),  flow  rate  (7.3  m/h) , 
FeCl,  suspension  (5  mg/I.,  pH=8.3). 

-1  -1 

Run  VII-2  O:  first  compartment  (0  "^   4.24  cm),  a=1.035  h   ,  B=0,06  h   , 

A:  second  compartment  (4.24  %  13.97  cm),  a=0.4133  h   , 

6=0.0286  h"''. 


,  -1 


third  compartment  (13.97  "^   44.45  cm),  a=0.089  h   ,  3=0, 
lumped  model  (0  "-^   44.45  cm),  a=0.469  h"  ,  B=D.076  h"  , 
fourth  compartment  (44.45  '^   60.96  cm). 
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Figure  V-7.   Fitting  the  present  model  to  Rimer's  data  (1968); 

bed  depth  (60.96  cm:  anthracite,  20.32  cm,  d  :  0.100  cm;  sand,  20.32  cm, 
d  =0.071  cm,  garnet  20.32  cm,  d  =0.059  cm),  flow  rate  (7.3  m/h) . 

s  8  _1  _1 

Run.VII-2;0:  first  compartment  (0  "^   13.97  cm),  a=0.778  h   ,  B-0.055  h   , 

A:  second  compartment  (13.97  "«   21.08  cm),  a=0.367  h   , 

S=0.126  h'""", 

;;    O:  third  compartment  (21.08  •^'  60.96  cm),  a=0.129  h"  ,  B=0, 

O:  lumped  model  (0  '^.  50.96  cm),  a-0.263  h"  ,  6=0.076  h'  . 
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Figure  V-8.   Fitting  the  present  model  to  Deb's  data  (1969): 

sand  size  (0.0647  cm),  bed  depth  (61  cm),  flow  rate  (A. 9  m/h) ,  Fuller' 
earth  suspension  (45  x 10"   vol%) . 

Run  I;  O:  first  compartment  (0  "^   12.2  cm),  a=0.261  h"  ,  8=0, 

A:  second  compartment  (12.2  ^  60  cm). 
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Figure  V-9.   Fitting  the  present  model  to  Ives'  data  (1961): 

sand  size  (0.0544  cm),  bed  depth  (20  cm) ,  flow  rate  (4.9  m/h) , 
chloreila  algae  suspension  (135  x  lO"^  vol%) . 

Run  IB;  O:  first  compartment  (0  "^   9.2  cm),  a=0.359  h   ,  B=0.102  h"  , 

A:  second  compartment  (9.2  cm  ^   20  cm). 


APPENDIX  A.   DERIVATION  OF  THE  DETERMINISTIC  GOVERNING  EQUATION,  EQ.  V-12,  IN 
THE  TEXT 

For  any  arbitrary  functions  f  and  g,  the  following  expression  holds 

(see,  e.g..  Van  Kampen,  1981): 

N-1  N         ^ 

r   g(n)   lEf(n)  -   J:  f(n)  TE   -"gCn)  (A-1) 

n=0  n=l 

For  the  case  where 

g(-l)  =  f(0)  =  g(N)  =  f(N+l)  =  0, 

Eq.  A-1  becomes 

N  N         _ 

Z      g(n)   lEfCn)  =  I    f(n)   E   g(n)  (A-2) 

n=0  n=0 

From  Eq.  V-11  In  the  text,  the  Master  equation  for  the  birth-death  process 
is  expressed  as 

dp  (c)  I 

—2 =   (E  -  Dp  p  (t)  +  (IE  •"-  1)X  p^(t)  (A-3) 

dt  n  n  n  n 

n  =  0.  1,  2,  ...  ,  n 


Multiplying  both  sides  of  Eq.  A-3  by  the  respective  n's  and  summing  all 
resultant  equations  from  n=0  to  n=n-,  we  obtain 

^  =   Z     n(]E-  Dp  p  (t)  +  Z     n(]E   -  DX  P  (t)         (A-4) 
•^"^      n=0  "  "      n=0 
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From  Che  property  of  the  one-step  operator  as  given  in  Eq.  A-2,  we  see  that 
Eq.  A-4  is  equivalent  to 

"O  ,         °0 

■^^21  >   Z  y  p  (t)(]E"-^-  l)n  +  Z     X  p  (t)(i:  -a)n         (A-5) 
<*'      n=0  "  °  n=0  "  " 

For  the  linear,  birth-death  process,  Che  intensities  of  transition,  X^  and 

p  ,  cake  Che  following  f onns ; 
n 

X   =   a(n.  -  n)  (A-6) 

n        0 

and 

n 
Substitution  of  Eqs.  A-6  and  A-7  into  Eq.  A-5  yields 

^=   -B<n>  +  c[n„  -  <n>]  (A-8) 

dt  lj 

This  expression  is  the  same  as  the  decerminiscic  governing  equacion,  Eq.  v-12, 
in  the  text,  if  <n>  is  replaced  with  n(c). 


APPENDIX  B.   EXPANSION  OF  ONE-STEP  OPERATOR,  EQS.  V-18  AND  V-19,  IN  THE  TEXT 
Since 

n  =  ngif  (t)  +"01; 

as  Indicated  by  Eq. 

one  is  equivalent  to  the  increase  of  ;;  by  n    .   Expanding  n  (c+n^   ,  t) 

in  the  Taylor's  series  form,  we  obtain 


-1/2  ..    .   ,_  .^  .„-l/2^  ,(;,0  +1   -l-^,(c,t)  +  ... 


(^-"o'''i4"o'^^--^^<=-^^      ^«-" 


This  in  conjunction  with  Eq.  V-16  in  the  text, 


1/2 
Tr(?,t)  =  n^   p(n,t), 


yields 


n-/%(n.l,t)  =  (1  .  n-/^  4  4  n-  -i^  .  ...)n-n(n.t) 


p(n+l,t)  =  (1  +  n-^/'  3I  +  I  n-^  ^  +  ...)p(n.t)  (B-2) 


lEp(n,t)=  p(n+l,  t) 

2 
.(l+n-l/^^^in-/^+...)p(n,t).  (B-3) 


3C 


Similarly, 


-1/2  _3  ^  1  ^-1  _1_  .      ,,  (B-5) 

"O  3?        2  "O      ,2 


APPENDIX  C.   TRANSFORMATION  METHOD  USED  IN  EQUATION  V-26  IN  THE  TEXT 
The  central-difference  approximation  gives 


^  f(,)  .  fC;  +'^0  -  HO  (c-i) 


_3i  f(^)  ,  f(c  +  AC)  -  2f(0  +  f(g  -  AQ  ((,-2) 

3C^  (AO^ 


Hence , 

=  ;^  [5:g(c)f(';  +  AC)  -  Eg(Of(c)]  Cc-3) 

From  the  property  of  the  one-step  operator  as  indicated  in  Eq.  A-2,  we 
obtain 


^[Ig(C)f(C  +  AC)  -  zg(c)f(c)] 


-^[Ef(c)8(c  -  AC)  -  J:f(c)g(c)] 


-z:f(c)  -57  g(c)  <c-.4) 

oC 


Eqs.  C-3  and  C-4  yield 


ZgCc)  ^  f(c)  =  -SfCc)  y|  g(C)  (C-5) 


Similarly,  we  have 


:>2  ,2 

z:g(c)  A  f(c)  =  mo  -^  g(c)  (c-6) 


is-ao 


APPENDIX  D.   JUSTIFICATION  OF  THE  INSTABILITY  OF  EQ.  V-3A  IN  THE  TEXT 
The  general  solution  of  Eq.  V-33  in  the  text  takes  the  form 

<i;>  =  c(l  +  ax)  (D-l) 

To  satisfy  the  initial  condition 

<C>  =  0   at  T  =  0, 

c  must  be  zero.   However,  it  is  postulated  that  since  the  solution  is  an 
approximation-,  a  slight  variation  possibly  will  exist  in  the  initial  value 
of  <(;>  so  that  c  may  take  a  small  value  other  than  zero.   Notice  that  for 
a  nonzero  c,  <C,>   grows  with  time  as  indicated  in  Eq.  D-l.   Thus,  a  nonzero 

c  implies  that  as  the  time  elapses,  <^>  will  increase  to  such  a  large  value 

1/2 
that  the  assumption  of  the  fluctuating  term  being  of  order  n,.   ,  as  given  in 

Eq.  V-15  in  the  text,  is  no  longer  valid.  It  is  worth  noting,  therefore,  that 
if  c  is  not  zero,  the  expansion  technique  employed  in  the  text  need  be  further 
justified  for  a  large  time. 

The  applicability  of  the  second-order,  pure  birth  model  developed  in  the 
text  is  examined  here.   For  the  filtration  process,  only  the  initial  operating 
period  is  of  interest;  the  operation  need  be  terminated  by  the  time  the  pre- 
ssure drop  buildup  becomes  too  high  or  the  filter  falls  to  achieve  the  required 


which  substantially  minimizes  the  fluctuations.   To  estimate  the  deviation 
involved  in  the  model  by  assuming  that  c=0.01,  let  us  consider  the  pressure 

drop  data  from  the  first  compartment  (0  ^   4.27  cm)  of  the  filter  of  Eliassen 

2 
(1935),   The  filter  bed  had  a  cross-sectional  area  of  1876.64  cm  .   The  poro- 
sity of  grains  packed  in  the  bed  has  been  determined  to  be  0.41,   The  diameter 
of  suspended  solids  is  0.00124  cm.   By  assuming  that  the  size  of  the  pores  is 
5  times  that  of  the  suspended  solids,  the  number  of  open  pores,  n,.,  is  then 
calculated  as 
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(4.27)(1876.64)(0.41)  _         11 

z -  fa. Do  X  10 

5Tr(0. 00124)  /6 

Obviously,  n  is  indeed  a  large  number  with  an  order  of  10  .   The  duration 

of  Run  6,  expressed  in  terms  of  dimensionless  time,  at,  is  12.3  as  given  in 

Figure  V-3  in  the  text.   The  mean  number  of  blocked  pores,  <n>,  consists  of 

the  macroscopic  and  fluctuating  terms  as  given  in  Eq.  V-28  in  the  text,  i.e., 

1/2 
<n>  =  ng<t)  +  Hg'  <5>  (D-2) 

Substituting  the  solution  for  ())  and  <;;>  as  given  in  Eqs.  V-32  and  V-36, 
respectively,  in  the  text  into  this  expression  gives 

"O"^       1/2        2 


mean  number  of  blocked  pores  at  the  end  of  filtration  as 

^"^  =  ^^'^VfSP^^  (0.01)  (6.58  X  10^^)1/^1  ^   12.3)2 


macroscopic  term  (first  term).   It  is  conceivable,  therefore,  that  even  toward 
the  end  of  filtration,  the  fluctuating  term  is  still  negligible  compared  with 
the  macroscopic  term.   Based  on  these  observations,  we  believe  that  the  second- 
order,  pure  birth  model  is  indeed  applicable. 
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APPENDIX  E.   DERIVATION  OF  EQUATION  V-44  IN  THE  TEXT 

The  differential  equation,  Eq.  A-8,  describes  the  linear,  birth-death 
process  as 

iiSt  =  -S<n>  +  c.[n„  -  <n>]  (E-D 

dt  U 

Solution  of  this  equation  subject  to  the  initial  condition 

<n>  =  0  at  t  =  0  (E-2) 


gives 


<n>  =  an 


-(a+e)t 
^  -  "^       )  (E-3) 

0    a  +  B 


Substitution  of  this  expression  Into  Eq .  V-40  in  the  text  yields 

,^\  g  +  B    1  (1-4) 

"^'^  =  "O^  ^   -(a+B)t'  ^ 

B  +  ae 

This  equation  in  conjunction  with  Eq.  V-42  in  the  text  results  in 

P 
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CHAPTER  VI   EXPERIMENTAL  STUDY  OF  DEEP  BED  FILTRATION: 
Stochastic  Analysis 

In  Chapters  IV  and  V,  stochastic  models,  namely,  the  linear,  birth- 
death  and  second-order,  pure  birth  processes,  coupled  with  the  Carman- 
Kozeny  equation,  have  been  employed  to  simulate  the  pressure  drop  and  its 
temporal  change  during  a  filtration  run. 

Performance  of  a  deep  bed  filter  depends  on  a  large  number  of  variables 
including  the  size  distribution  of  suspended  solids,  suspension  concentration, 
filtration  rate,  and  bed  characteristics.  The  latter  two  variables  have  been 
widely  investigated  by  numerous  researchers  (Rimer,  1968;  Deb,  1969;  Tchobano- 
glous,  1970;  Tchobanoglous  and  Eliassen,  1970;  Huang,  1972).  The  primary 
concern  of  the  work  presented  in  this  chapter  is  to  Investigate  experimentally 
the  effect  of  the  size  distribution  and  concentration  of  suspended  solids  on 
the  pressure  drop  dynamics  within  the  filter  bed;  the  results  are  treated 
stochastically. 

EXPERIMENTAL  FACILITIES  AND  METHODS 

A  schematic  of  the  experimental  set-up  is  shown  in  Figure  Vl-1.   The 

set-up  comprised  an  upflow  deep  bed  filter,  a  suspension  tank,  a  water  tank, 

two  pumps,  a  differential  pressure  transducer,  a  transducer  indicator,  and  a 

strip  chart  recorder.  The  plexiglass  filter  had  an  Inside  diameter  of  0.132  m 

and  a  length  of  0.406  m  between  the  bottom  distributor  plate  and  the  outlet. 

The  filtering  medium  consisted  of  -  #20  to  +  #30  (20x30)  mesh  silica  sand  with 

3 
an  average  diameter  of  approximately  0.710  mm  and  a  density  of  2600  kg/m  . 

The  bed  of  sand,  5  cm  in  depth,  was  supported  between  a  bottom  distributor 

plate  and  a  top  movable  porous  plate.  These  plates  retained  the  sand 

particles  but  allowed  the  suspension  to  flow  through. 

The  suspension  was  composed  of  coal  particles  dispersed  in  water. 
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Five  different  size  ranges  of  coal  particles,  which  were  carefully  sieved 
at  least  three  times,  were  used  in  this  work;  these  size  ranges  are  outlined 
in  Table  VI-1.  For  each  specific  size  range  of  coal  particles,  the  suspension 
concentration  ranged  approximately  from  0.001  to  0,015  weight  percent  (wtt,) . 
The  coal  particles  were  maintained  in  water  in  the  suspension  tank  by  agitating 
it  with  a  1/4  HP  twin  propeller  mixer. 

The  porosity  of  the  bed  was  dependent  on  the  packing  conditions;  it 
ranged  from  0.43  under  the  least  dense  static  condition  to  0.35.  Prior  to 
a  filtration  run,  the  sand  was  allowed  to  settle  and  form  a  bed  with  a  porosity 
of  around  0.38  with  the  top  porous  plate  placed  and  maintained  at  a  position 
immediately  above  the  sand.  The  influent  concentration,  pH  and  temperature  of 
the  suspension  were  measured  before  each  filtration  run.   Suspension  from  the 
suspension  tank  was  introduced  at  a  constant  rate  of  48.8  m/h  to  the  bottom 
of  the  filter  by  a  HYPRO  piston  pximp.  The  suspension  was  then  allowed  to  flow 
upward  through  the  filter. 

The  pressure  drop  increase  in  the  filter  bed  was  monitored  with  a  vari- 
able reluctance  pressure  transducer  which  was  connected  to  pressure  taps  posi- 
tioned below  the  bottom  distributor  plate  and  near  the  outlet  of  the  filter, 
as  shown  in  Figure  VI-1.  A  strip  chart  recorder  connected  to  the  pressure 
transducer  continuously  recorded  the  change  in  the  pressure  drop  during  a  fil- 
tration run.  The  pressure  drop  attributable  to  the  column  wall  and  porous 
plates  was  measured  by  passing  water  through  the  empty  column  at  the  same  rate 
of  48.8  m/h;  it  was  determined  to  be  1.73  cm  H„0.  This  amount  was  subtracted 
from  the  recorded  pressure  drop  to  obtain  the  pressure  drop  through  the  filter 
bed.  Samples  of  the  filtrate  were  collected  intermittently  from  the  effluent 
flowing  through  the  outlet  without  disturbing  the  performance  of  the  filter. 
These  samples  were  then  filtered  in  a  200  mSL  pressure  filter  to  determine  the 
solids  content. 


VI-3 
The  filter  was  cleaned  by  fluldizing  the  entire  bed  with  water  to  remove 
the  deposited  solids.  During  this  period,  the  upper  porous  plate  was  raised 
to  a  position  well  above  the  fluidized  sand  particles.   The  coal  particles 
were  easily  flushed  from  the  filter  in  this  manner.  With  the  bed  regenerated, 
the  filtration  cycle  could  be  reinitiated. 

RESULTS  AND  DISCUSSION 

Five  sets  of  filtration  runs,  namely,  runs  of  series    A,  B,  C,  D,  and 
E,  were  conducted  with  suspended  solids  of  various  sizes,  size  distributions, 
and  concentrations.   For  each  run,  the  pressure  drop  and  the  effluent  concentra- 
tion were  measured  as  functions  of  filtration  time. 

Modeling  of  Pressure  Drop  Dynamics 

Stochastic  models  which  simulate  the  pressure  drop  dynamics  within  the 
filter  bed  under  the  constant  flow  condition  are  categorized  into  two  classes, 
the  linear,  birth-death  and  second-order,  pure  birth  processes,  as  presented 
In  Chapters  IV  and  V,  respectively. 

Models.   The  linear,  birth-death  model  assumes  that  the  blockage  Intensi- 
ty. 1  ,  and  the  scouring  intensity,  u  ,  are  linearly  proportional  to  the  number 
^  '  "'n '  '^  n 

of  open  pores  and  the  number  of  blocked  pores  In  the  bed,  respectively  (Lltwini- 
szyn,  1968a;  1968b;  1969).  This  assumption  of  linearity  gives  rise  to 

X  =  a(n„  -  n);  n  =  0,  1,  2,  . . . ,  n         (VI-1) 

n      u  ^ 

/i  =^n;  n  =  0,  1,  2,  ...,  Ug         (VI-2) 

where  n„  is  the  total  number  of  open  pores  susceptible  to  blockage  In  the 
bed,  a   Is  the  blockage  constant,  and  p   is  the  scouring  constant.  As  presented 
In  Chapter  IV,  the  mean  number  of  blocked  pores  at  the  moment  t,  <n> ,  can  be 
solved  based  on  Che  assumption  given  in  Eqs .  VI-1  and  VI-2.   Relating  the 
mean  niimber  of  blocked  pores  to  the  linear  or  interstitial  velocity  of  the 
flow  and  incorporating  the  resultant  expression  into  the  Carman-Kozeny 
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equation  give 

a  B  +  ae 

P 

Under  Che  simplifying  assumption  that  the  increase  in  the  pressure  drop  can 
be  attributed  to  that  of  the  linear  velocity  of  the  flow,  we  have 


(.AZii),  .(_^)(.^iJ__,  (vi-A) 


-^HOl,  .JL±1 (VI-5) 


where  r-AP(t)]  represents  the  pressure  drop  through  the  bed  at  moment,  t, 
and  f-^APJ  that  at  the  onset  of  filtration.   Notice  that  if  the  scouring  is 
negligible,  i.e., 


^  =  0 
the  linear,  birth -death  process  reduces  to  the  linear,  pure  birth  process 
(Litwiniszyn,  1963;  Hsu  and  Fan,  1984).   The  resultant  equation  is 

The  second-order,  pure  birth  model  assumes  that  X      is  proportional  to 
the  square  of  the  number  of  open  pores  in  the  bed,  and  fx     is  negligible. 
This  assumption  gives  rise  to 

^  =  °  ■        .  .  \'      - 

As  presented  in  Chapter  V,  the  Master  equation  governing  the  probability 
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balance  can  be  solved  approximately  by  means  of  the  so-called  expansion 
technique.   This  technique  yields  the  mean  number  of  blocked  pores  and 
subsequently   results  in  the  following  equation; 

1-  AP(t)l  ^  1  .  .,  (VI-9) 


k  =  an  <VI-10) 

k   unQ 

Comparison  of  pressure  drop  data  with  models.      The   present  models 
consider  the  entire  bed  as   a  single  spatially-lumped  compartment.      For  the 
filter  bed  with  a  depth  of  merely  5  cm,   it   is   conceivable   that  such   lumped 
models  are  satisfactory. 

In   Runs  A-1    through  A-9,    the    suspension   filtered  was   composed   of 
-  #270    to  +  #400    (270x400)    mesh   coal   particles  with  an  average   diameter  of 
approximately  0.044  mm.      The   pressure  drop  was   measured  as   a   function  of 
filtration    time    for  each   run;    the    results    of   Run  A-1,  A-4,    and  A-8   are   shown 
in  Figure  VI-2.      This    figure   shows   that   the  change   in   the   pressure  drop 
through   the   bed   can  be   described    reasonably  well   by   the-  linear,    birth- 
death  model    (Eq.    VI-5) . 

The   suspension,    used  in  Runs   B-1    through    B-9,    comprised    -  #250    to 

+  #270    (250x270)    mesh   coal    particles  with   an  average  diameter  of  approximately 
0.058  mm.      The    pressure   drop  data    from  Runs    B-2,    B-5,    and    B-9  are    plotted 
against    the    filtration    time    in   Figure   VI-3.      This    figure  demonstrates    that 
the   change    in    the    pressure  drop    through    the   bed   also   exhibits    the    linear, 
birth-death  behavior. 

The    suspension,    filtered    in   Runs   C-1    through  C-6,    consisted   of    -  #140 
to  +  #200    (140x200)    mesh   coal   particles   with   an  average   diameter  of  approxi- 
mately  0.090  mm.      The    pressure   drop  data    from  Runs   C-1,   C-4,    and  C-6   are 


■   ,vs^:,  •  VI-6 

plotted  in  Figure  VI-4,   This  figure  indicates  that  the  change  in  the 
pressure  drop  through  the  bed  obeys  essentially  the  second-order,  pure 
birth  process  (Eq.  VI-9) . 

In  Runs  D-1  through  D-6,  the  suspension  filtered  was  composed  of 
-  #100  to  +  #140  (100x140)  mesh  coal  particles  with  an  average  diameter  of 
approximately  0.127  mm.   The  pressure  drop  was  measured  as  a  function  of 
filtration  time  for  each  run;  the  results  of  Runs  D-1,  D-4,  and  D-6  are 
shown  in  Figure  VI-5.  This  figure  shows  that  the  change  in  the  pressure 
drop  during  the  initial  period  of  a  filtration  run  reflects  the  linear, 
pure  birth  behavior  (Eq.  VI-6).  However,  it  also  shows  a  departure  from 
the  predicted  behavior  in  the  later  stage  of  the  filtration  run.   Note 
that  the  diameter  ratio  of  coal  particles  to  grains  is  up  to  0.18,  indicating 
that  a  surface  mat  of  deposit  will  form  on  the  bottom  of  the  filter  during 
the  later  stage  of  a  filtration  run.  This  implies  that  the  deep  bed 
filtration  will  tend  to  resemble  the  cake  filtration;  the  entire  bed  simply 
serves  as  the  meditmi  for  the  cake  formation  of  coal  particles  during  this 
period.  Therefore,  the  pressure  drop  increase  will  certainly  be  greater 
that  predicted  for  the  deep  bed  filtration.   Naturally,  this  phenomenon  of 
cake  formation  is  undesirable  and  should  be  avoided  in  the  design  of  the 
deep  bed  filter. 

In  Runs  E-1  through  E-13 ,  the  suspension  filtered  consisted  of  -  #100 
to  +  #270  (100x270)  mesh  coal  particles  with  an  average  diameter  of  approxi- 
mately 0.101  ram;  that  is,  a  wider  size  range  of  coal  particles  was  employed 
in  the  filtration  runs  to  observe  the  combined  effect  of  all  possible 
mechanisms  predicted  by  the  present  models.   It  appears  that  the  combined 
mechanism  of  the  linear,  pure  birth  process  and  cake  filtration  predominates 
over  either  the  linear,  birth-death  or  second-order,  pure  birth  mechanism 
as  shown  in  Figure  VI-6,  Notice  that  the  pressure  drop  Increase  in  the 
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filter  bed  tends  to  be  larger  than  that  in  any  of  the  previous  runs  for  the 
same  suspension  concentration.   Similarly,  the  capacity  of  the  filter  with 
unisized  grains  is  greater  than  that  with  grains  varying  widely  in  size; 
that  is,  the  pressure  drop  increase  in  the  former  is  smaller  than  that  in 
the  latter  (Streander,  1940).   These  results  imply  that  the  size  uniformity 
of  grains  or  suspended  solids  may  reduce  the  tendency  of  suspended  solids  to 
clog  the  bed,  thus  reducing  the  rate  of  pressure  drop  increase. 

Analysis  of  Effluent  Concentration  Data 

The  effluent  concentrations  of  suspended  solids  for  Runs  A-1,  A-2,  A-5, 
and  A-8  are  plotted  as  functions  of  filtration  time  In  Figure  VI-7.   Note 
that  the  effluent  concentrations  increase  with  filtration  time.   This  is 
indeed  indicative  of  the  existence  of  the  scouring  phenomenon  within  the  bed; 
recall  that  the  scouring  rate  increases  with  an  increase  in  the  deposit 
as  indicated  in  Eq.  VI-2.   This  monotonically  increasing  trend 
is  also  observed  in  the  effluent  concentration  data  for  Run  B-1  as  shown  in 
Figure  VI-8.  Therefore,  the  use  of  the  linear,  birth-death  model  is  more 
than  adequate  for  representing  the  pressure  drop  dynamics  in  Runs  A  and  B. 

For  Run  C-1,  the  effluent  concentration  of  suspended  solids  was  measured 
as  a  function  of  filtration  time  as  shown  in  Figure  VI-9.  Notice  that  the 
data  do  not  show  the  monotonically  increasing  trend;  they  seem  to  be  inde- 
pendent of  filtration  time.  Thus,  the  second-order,  pure  birth  model,  for 
which  the  scouring  mechanism  is  assumed  to  be  negligible,  need  be  employed 
to  depict  the  pressure  drop  dynamics  in  Runs  C-1  through  C-6. 

The  effluent  or  filtrate  obtained  from  each  of  Runs  D-1  through  D-6 
was  essentially  clean  water.   This  implies  that  during  a  filtration  run 
the  rapid  clogging  of  suspended  solids  around  the  entrance  section  of  the 
bed  will  be  significant.   It  is  plausible,  therefore,  that  the  cake  filtration 
dominated  in  these  filtration  runs. 
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Effect  of  Suspension  Concentration 

The  pressure  drop  dynamics  in  Runs  A-1  through  A-9  and  in  Runs  B-1 
through  B-9  can  be  depicted  by  the  linear,  birth-death  model  with  two 
parameters,  namely,  the  blockage  constant,  a,  and  the  scouring  constant,  j3. 
The  suspension  concentration  ranged  from  0.0147  wt%  for  Run  A-1  to  0.0015  wt7„ 
for  Run  A-9,  and  it  ranged  from  0.0084  wt%  for  Run  B-1  to  0.0013  wt7o  for 
Run  B-9.  The  dependence  of  a  and  /3  on  the  suspension  concentration  for 
the  runs  of  series  A  is  illustrated  in  Figure  VI-10  and  that  for  the  runs 
of  series  B  is  shown  in  Figure  VI-11.   Both  figures  indicate  that  a  and  |8 
increase  with  an  increase  in  the  suspension  concentration.  Notice  that  /3 
drops  sharply  to  approach  zero.  This  implies  that  the  blockage  mechanism 
will  predominate  over  the  scouring  mechanism  while  the  suspension  becomes 
sufficiently  dillute,  thus  giving  rise  to  the  linear,  pure  birth  process. 

The  second-order,  pure  birth  model  is  employed  to  simulate  the  pressure 
drop  characteristics  in  Runs  C-1  through  C-5.  The  suspension  concentration 
employed  ranged  from  0.0101  wt7.  for  Run  C-1  to  0.0018  wt%  for  Run  C-6.  The 
dependence  of  the  constant,  k,  defined  in  Eq.  VI-10,  on  the  suspension 
concentration  is  shown  in  Figure  VI-12.   It  appears  that  k  also  increases 
with  an  increase  in  the  suspension  concentration. 

The  linear,  pure  birth  model  can  be  employed  to  describe  the  change 
in  the  pressure  drop  during  the  earlier  period  of  a  filtration  run  in  Runs 
D-1  through  D-6  and  in  Runs  E-1  through  E-13.  The  suspension  concentration 
ranged  from  0.0104  wt7.  for  Run  D-1  to  0.0036  wt7o  for  Run  D-6,  and  it  ranged 
from  0.0122  wtT,  for  Run  E-1  to  0.0031  wt7o  for  Run  E-13.  The  dependence  of 
the  blockage  constant,  a,  on  the  suspension  concentration  for  the  runs  in 
series  D  is  depicted  in  Figure  VI -13,  and  that  for  the  runs  of  series  E 
is  shown  in  Figure  VI-14.   Both  figures  show  that  a  increases  with  an 
increase  in  the  suspension  concentration. 
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Effect  of  Size  and  Size  Distribution  of  Suspended  Solids 

Figure  VI-15  illustrates  the  effect  of  the  size  distribution  of 
suspended  solids  on  the  dynamic  behavior  of  pressure  drop  through  the 
bed  for  the  same  suspension  concentration  of  0.0073  wt%.  As  the  particle 
size  increases,  the  pressure  drop  dynamics  undergo  transition  from  the 
linear,  birth-death  behavior,  through  the  second-order,  pure  birth  behavior, 
to  the  linear,  pure  birth  behavior;  they  eventually  exhibit  the  characteris- 
tics of  the  cake  filtration.  Notice  that  the  smaller  the  size  of  suspended 
particles,  the  easier  they  escape  from  the  bed.   It  can  be  postulated  that 
-relatively  speaking,  the  small  particles,  once  deposited  in  the  bed,  will 
tend  to  be  easily  scoured  by  the  forces  acting  on  them.  On  the  contrary, 
the  large  suspended  particles  tend  to  be  trapped  within  the  bed.   This 
implies  that  as  the  size  of  suspended  particles  increases  to  a  certain 
degree,  the  particles  will  be  unable  to  penetrate  into  the  bed  to  any  great 
extent  before  being  captured  by  the  grains.  And  thus,  a  "cake"  of  suspended 
particles  will  begin  to  build  up,  and  a  sharp  increase  in  the  pressure  drop 
will  be  observed  as  can  be  seen  in  Figures  VI-5  and  VI-6. 
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Figure  VI-2.   Increase  in  the  pressure  drop  as  a  function  of  filtration  time  : 

"^,  S  =  0.1277  min''",  pH  =  9.1,  T  =  20.6°C), 


Runs;  A-1   (a  =  0.0324  min 
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0.0374  min"  ,  pH  =  9.1,  T  =  20.7  C). 
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Figure  VI-3.      Increase    In    the   pressure   drop  as   a    function   of    filtration    time    : 

Runs;    U-2      (a   =   0.0593  min''',    B   =   0.1256  min"    ,    pH   =   9.0,   T   =   20.6°C), 


B-5      (a   =  0.0409  min"',    B   =  0.0985  min"'' ,    pH   =   9.1,   T  =  20.8*C), 
B-9      (a   =  0.0159  min"-',    B   =  0.0102  min"',    pH   =  9.1,   T  =   20.9°C). 
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Figure  Vl-4.      Increase    in    Che    pressure   drop  as  a  function   of    filtration    time 

(k   =  0.0575   min"    ,  pli  =  8.5,    T   =    21.0°C), 

(k    =    0.0321    min"    ,  pll  =  9.0,    T    =    19.6  'c), 

(k   =   0.0174   min"',  pll  =  9.1,    T   =    18.7  'c). 
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Figure  VI-5.     Increase   in    Che  pressure   drop  as    a    function   of    filtration   time 

Runs;    D-1      (a    =  0.040  min"    ,    pH   =  8 . 7 ;    T   =   16.5 'c), 

D-4       P    =  0.027   min''"',    pH    =   8.9;    T    =    16.3    C)  , 

D-6      ("    =  0.0162  min'-'-,   pH   =  8.8;    T  =  16.3  C) . 
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Figure  Vl-6.   Increase  in  the  pressure  drop  as  a  function  of  filtration  time  : 

Runs;  E-1  (a   =  O.lOOO  min''' ,  pH  =  8.3,  T  =  21.9  °C), 

E-6   (O  =  0.0579  min'  ,  pH  =  8.6,  T  =  21.7  °C), 

E-10  (o  =  0.0357  mln'''^,  pH  =  9.2,  T  =  22.3  °C). 


VI-17 


to 


o 

X 


c 
o 


c 

CD 
O 

a. 
o 
O 

c 

QJ 

Z3 


LJ 


15 

14 

13- 

12- 

11 

10- 

9- 

8 

7 

6 

5 

4 

3 
2 

1 
0 


O 


A    -   1 

o    o    o    o 


o 


A      A 


O 


O    o     o 


□       °       ° 
J [ L 


O 


C      =  0.0147  wt7. 

L 

A  C.    =  0.0133  wt7. 


A   -  5 

^       O  C.    =  0.0061  wt7. 


A    -  8 

D        □  C     =  0.0021  wt7. 

1 

I       1       I I \ L 


0      1      2      3     4      5      6      7      8      9     10    11     12    13    14    15 

Filtration  Time(min) 

Figure  VI-7.   Effluent  concentration  of  suspended  particles  as  a  function 

of  filtration  time  :  Runs  A. 


VI-18 


O 

X 


o 


c 

O) 

o 

o 
(J> 

c 

33 


LJ 


3       4        5        6       7        8        9 

Filtration  Time  (min) 


10     11      12 


Figure  VI-8. 


Effluent   concentration   of   suspended    particles   as   a    function 
of   filtration  time    :    suspension  concentration,    0.0084  wt%; 
Runs    B-1. 
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Figure  VI -9. 


Effluent  concentration  of  suspended  particles  as  a  function 
of  filtration  time  :  suspension  concentration,  0.0101  wtX; 
Runs  Crl. 
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Figure  VI-10.   Dependence  of  the  blockage  constant,  a,    and  the  scouring  constant, 
g,  on  the  suspension  concentration  :  Runs  A -1  through  A-9. 
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Figure  VI-11.  Dependence  of  the  blockage  constant,  a,  and  the  scouring  constant, 
R,  on  the  suspension  concentration  :  Runs  B-1  through  B-9. 
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Figure  VI-12.   Dependence  of  the  constant,  k ,  on  the  suspension  concentration 
Runs  C-1  through  C-6. 
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Figure  VI-13.   Dependence  of  the  blockage  constant,  a,  on  the  suspension 
concentration  :  Runs  D-1  through  D-6. 
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Figure  VI-14.   Dependence  of  the  blockage  constant,  a   ,  on  the  suspension 
concentration  :  Runs  E-1  through  E-13. 
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Figure  VI-15.  Increase  in  chc  pressure  drop  as  a  function  of  filtracion 
time  with  Che  size  of  coal  particles  as  the  parameter  : 
suspension   concentration  ,    0  .  0073    vt'A. 


CHAPTER  VII  AN  EXPERIMENTAL  STUDY  OF  SEMIPLUIDIZED- BED  FILTERS: 

Filtration  of  Suspended  Solids  with  Narrow  Size  Ranges 

Deep  bed  filters,  also  known  as  granular  filters,  are  used  extensively 
in  purification  of  potable  water  and  in  treatment  of  waste  water.   The  oper- 
ation of  deep  bed  filters  is  inefficient,  however;  It  is  typically  character- 
ized by  the  shallow  penetration  of  suspended  solids.  As  discussed  in  Chapter 
IV,  Ives  (1961)  measured  in  detail  the  specific  deposit  as  a  function  of  time 
and  bed  position  by  means  of  the  radioactive  algae  technique  and  concluded 
that  the  penetration  of  suspended  solids  through  his  deep  bed  filter  was  indeed 
shallow.   The  sectional  analysis  of  the  pressure  drop  dynamics  within  the  deep 
bed  filter  presented  in  Chapter  V  indicates  that  the  pressure  drop  through  the 
last  portion  of  the  filter  bed  usually  remains  invariant  with  time  and  this 
unutilized  portion  can  be  as  much  as  80%  of  the  entire  bed  as  illustrated  in 
Figure  V-8  (Deb,  1969).   Furthermore,  a  layer  of  filtered  solids  will  usually 
collect  on  the  surface  and  within  the  first  few  layers  of  the  grains;  that  is, 
a  "cake"  of  suspended  solids  will  be  formed  around  the  entrance  section  of  the 
bed,  thus  giving  rise  to  an  abrupt  increase  in  the  pressure  drop.   This  pheno- 
menon of  "cake"  formation  has  been  discussed  in  Chapter  VI  (see  Figures  VI-11 
and  VI-13).   The  shallow  penetration  and  "cake"  formation  lead  to  Inefficient 
and  short  operation  of  the  deep  bed  filters,  thus  requiring  frequent  backwashing 
(see,  e.g..  Diaper  and  Ives,  1965;  Rimer,  1968;  Huang,  1972). 

To  enhance  the  penetration  of  suspended  solids  into  the  filter,  the  size- 
graded  sand  filter  was  introduced  to  replace  the  unisized  granular  filter.   The 
down-flow  filter  comprises  sand  of  decreasing  size  in  the  direction  of  flow. 
The  clean  filter  performs  well;  however,  during  the  backwash  cycle,  the  larger 
grains  tend  to  settle  to  the  bottom  of  the  bed  while  the  smaller  ones  tend  to 
rise  to  the  top,  resulting  in  a  reverse  stratification  in  the  filter.   This 
eventually  causes  the  filter  performance  to  deteriorate  (Diaper  and  Ives,  1965). 
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Reverse  stratification  of  collector  grains  can  be  prevented  by  employing 
a  multimedia  filter.   The  multimedia  filter  is  composed  of  layers  of  filter- 
ing media;  the  particle  size  of  each  succeeding  layer  decreases  similar  to 
the  size  graded  filter,  while  the  density  of  the  media  increases  with  each 
succeeding  layer.   The  larger  but  lighter  particles  remain  at  the  top  of  the 
bed  and  the  smaller  but  heavier  particles  remain  at  the  bottom  during  and 
after  the  backwash  cycle.   A  typical  multimedia  filter  consists  of  anthracite, 
sand,  and  garnet  particles  (Diaper  and  Ives,  1965;  Rimer,  1968).   Huang  (1972) 
showed  that  filtration  runs  with  a  dual  media  filter  of  anthracite  and  sand 
were  60  to  160%  longer  than  those  obtained  with  a  unisized  sand  filter  were 
attainable.   Similarly,  Rimer  (1968)  showed  that  the  average  net  yield  was 
approximately  50%  higher  for  water  processed  in  the  anthracite-sand-garnet 
filter  than  in  the  conventional  sand  filter.   The  results  of  these  studies 
indicate  that  the  performance  of  the  sand  filter  is  improved  when  a  greater 
portion  of  the  bed  is  utilized  through  the  deeper  penetration  of  solids. 

The  semifluidized-bed  filter  proposed  by  Wen  and  Fan  (1979)  has  been 
employed  to  overcome  the  problem  of  the  rapid  clogging  around  the  entrance 
section  of  the  bed  or  the  phenomenon  of  "cake"  formation,  which  is  usually 
experienced  in  the  deep  bed  filtration.   Unlike  the  multimedia  filters,  only 
a  single  medium,  sand,  is  necessary.   Initially,  the  bed  will  be  in  a  completely 
fixed  configuration  as  the  fluid  stream  containing  suspended  solids  enters  the 
bed.   The  suspended  solids,  upon  entering  the  portion  of  the  bed  around  the 
entrance  section,  will  be  unable  to  penetrate  the  bed  to  any  great  extent 
before  being  stopped  by  the  grains  of  sand.   A  "cake"  of  suspended  solids  will 
begin  to  be  built  up  in  this  zone.   The  "cake"  formation  will  produce  a  sharp 
increase  in  the  pressure  drop.   This  rise  in  the  pressure  drop  will  activate 
a  driving  mechanism  which  will  raise  the  upper  restraining  plate.   The  mechani-: 
cal  movement  of  the  restraining  plate  in  conjunction  with  the  forces  being 
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exerted  by  the  fluid  stream  will  cause  disruption  of  the  formed  "cake"  of 
suspended  solids  into  small  fragments.   Cake  disruption  will  also  Ua  aided  by 
the  fluidized  particles.   Raising  of  the  restraining  plate  and  the  disruption 
of  the  formed  "cake"  will  continue  until  the  breakthrough  of  the  suspended 
solids  occurs,  which  means  the  ever  decreasing  fixed  bed  can  no  longer  achieve 
the  required  removal  of  the  suspended  solids  from  the  suspension. 

Another  type  of  the  seraifluidized-bed  filter  has  been  proposed  by  Hsu 
and  Fan  (1981;  1982)  and  Fan  and  Hsu  (1984)  to  overcome  the  problems  of  shallow 
solids  penetration  and  rapid  clogging  simultaneously.   In  operating  this  filter, 
the  filtration  cycle  starts  with  the  restraining  plate  in  a  position  forming 
the  desired  proportion  of  fixed  to  fluidized  beds.   During  the  filtration  run, 
the  fixed  section  functions  as  a  deep  bed  filter  in  the  capture  of  suspended 
solids.   The  fluidized  section,  on  the  other  hand,  provides  the  necessary 
turbulence  in  scouring  the  deposited  solids,  thus  delaying  the  pressure  drop 
increase. 

The  present  work  aimed  to  continue  the  investigation  on  the  semif luidized- 
bed  filter  proposed  by  Hsu  and  Fan  (1981;  1982)  and  Fan  and  Hsu  (1984).  However, 
instead  of  using  suspended  solids  varying  widely  in  size,  such  as  sub  #50  mesh 
coal  particles  used  in  their  work,  coal  particles  with  several  narrow  size 
ranges  were  employed;  these  include  -  #270  to  +  /MOO  (270x400)  mesh,  -  #250  to 
+  //270  (250x270)  mesh,  -  #140  to  +  #200  (140x200)  mesh,  and  -  #100  to  +  #140 
(100x140)  mesh  coal  particles.  The  primary  interest  was  to  investigate  the 
effects  of  both  the  size  of  suspended  solids  and  the  "cake"  formation  on  the 
performance  of  semif luidized-bed  filters, 

EXPERIMENTAL  FACILITIES  AND  METHODS 

A  schematic  of  the  experimental  set-up  is  shown  in  Figure  VIT-1.   The 
set-up  consisted  of  a  semif luidized-bed  filter,  a  suspension  tank,  a  water 
tank,  two  pumps,  a  differential  pressure  transducer,  a  transducer  indicator. 
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and  a  strip  chart  recorder.   The  semlfluldized-bed  filter  was  constructed  from 
a  plexiglass  column  with  an  inside  diameter  of  0.132  m  and  a  length  of  0.406  m 
between  the  bottom  distributor  plate  and  the  outlet.  The  filtering  medium  was 

composed  of  -  i^2Q   to  +  #30  (20x30)  mesh  silica  sand  with  an  average  diameter  of 

3 
approximately  0.710  mm  and  a  density  of  2600  kg/m  . 

The  minimum  fluidization  velocity  was  determined  experimentally  to  be 
13.1  m/h.   During  a  filtration  run,  suspension  from  the  suspension  tank  was 
introduced  at  a  constant  rate  of  48.8  m/h  to  the  bottom  of  the  filter  by  a  HYPRO 
piston  pump.   For  the  operation  of  the  deep  bed  filtration,  the  sand,  20  cm  in 
depth  (between  0.149  m  and  0.349  m  of  the  column)  with  the  porosity  of  around 
0.40,  were  supported  between  a  bottom  distributor  plate  and  an  upper  movable 
restraining  plate.   These  porous  plates  retained  the  sand  but  permitted  the 
suspension  to  pass.   This  same  amount  of  sand  also  served  as  the  filtering  medium 
for  the  operation  of  semif luidized-bed  filters. 

The  suspension  comprised  coal  particles  dispersed  in  water.   Coal  particles 
of  four  different  size  ranges,  which  were  carefully  sieved  at  least  three  times, 
were  employed;  these  size  ranges  are  listed  in  Table  VII-1.   During  a  filtration 
run,  the  concentration  of  suspension  was  maintained  constant  in  the  suspension 
tank  with  agitation  provided  by  a  1/4  HP  propeller  mixer. 

By  positioning  the  restraining  plate  at  various  heights,  different  ratios 
of  the  fixed-bed  to  f luidized-bed  heights  could  be  obtained.  Prior  to  starting 
a  filtration  run  under  the  constant  influent  concentration  and  flow  condition, 
the  upper  restraining  plate  was  moved  to  a  predetermined  position,  which  would 
result  in  a  desired  degree  of  fluidization  with  the  introduction  of  suspension. 
During  a  filtration  run,  the  upper  restraining  plate  was  maintained  at  this  posi- 
tion. For  each  suspension  containing  coal  particles  of  specific  size  range, 
filtration  with  the  semifluidized  bed  was  performed  with  75%  or  50%  fluidization, 

i.e.,  h  /h  =  25%  or  50%,  respectively,  in  which  h   represents  the  initial 
'   pa  >    r       J  >  pa 
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fixed-bed  height  and  h  the  entire  bed  height.  Filtration  runs  with  a  complete- 
ly fixed  bed,  namely,  a  20  cm  conventional  deep  bed  filter,  were  also  conducted 
under  the  same  influent  concentration  and  flow  condition  to  provide  a  basis  for 
comparison  and  evaluation. 

The  pressure  drop  increase  through  the  filter  bed  was  monitored  with  a 
variable  reluctance  pressure  transducer  which  was  connected  to  pressure  taps 
positioned  below  the  bottom  distributor  plate  and  near  the  outlet  of  the  filter, 
as  shown  in  Figure  VII-1.  A  strip  chart  recorder  connected  to  the  pressure 
transducer  indicator  continuously  recorded  the  change  in  the  pressure  drop 
during  a  filtration  run.  The  pressure  drop  attributable  to  the  column  wall 
and  porous  plates  was  measured  by  passing  water  through  the  empty  column  at  a 
constant  rate  of  48.8  m/h;  it  was  determined  to  be  1.73  cm  H„0.   This  amount 
was  subtracted  from  the  recorded  pressure  drop  to  obtain  the  pressure  drop 
through  the  filter  bed.   Samples  of  the  filtrate  were  collected  intermittently 
from  the  effluent  flowing  through  the  outlet  without  disturbing  the  performance 
of  the  filter.   These  samples  were  then  filtered  in  a  20  ml   pressure  filter 
to  determine  the  solids  content. 

The  exhausted  filter  was  cleaned  by  fluidizing  the  entire  bed  with  water 
to  remove  the  deposited  solids.   During  this  period,  the  upper  restraining  plate 
was  raised  to  a  position  well  above  the  fluidized  sand  as  seen  in  Figure  VII-2, 
The  coal  particles  were  easily  flushed  from  the  filter.  With  the  bed  regener- 
ated, the  filtration  cycle  could  be  re-activated, 

RESULTS  AND  DISCUSSION 

The  results  obtained  with  the  semifluidized-bed  filters,  namely,  the 
results  from  runs  of  series  F,  G,  H,  and  I,  are  presented  and  compared  to  those 
obtained  with  the  equivalent  deep  bed  filters  in  Figure  VII-3  through  VII-14. 
Also  presented  are  two  sets  of  photographs  taken  during  runs  with  two  semifluid- 
ized-bed filters.   One  set,  including  Figures  VII-9  and  VII-11,  is  for  the 
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filter  with  h  /h  =  25%  (Run  1-2),  and  the  other,  including  Figures  VII-12 

and  VII-14,  for  the  filter  with  h   /h  =  50%  (Run  1-3). 

pa 

Filtration  of  Relatively  Small  Particles 

In  Runs  F-1  through  F-3,  the  suspension  was  composed  of  -  //270  to  +  #400 

(270x400)  mesh  coal  particles  with  an  average  diameter  of  approximately  0.044  mm. 

The  influent  concentration  was  maintained  at  0.0227  wt%  for  each  filtration  run. 

The  increase  in  the  pressure  drop  across  the  entire  bed  was  the  largest  for  the 

deep  bed  filter  (h  /h  =  100%:  Run  F-1),  and  the  smallest  for  the  semif luidized- 
pa        '         ' 

bed  filter  with  75%  fluidization  (h  /h  =  25%;  Run  F-2)  as  seen  in  Figure  VII-3. 

pa 

The  effluent  concentration  for  all  runs  are  plotted  against  the  filtration 
time  in  Figure  VII-4.  The  steadily  increasing  effluent  concentrations  are  app- 
arently indicative  of  the  significant  scouring  effect  in  these  runs;  the  only 
exception  is  Run  F-2,  in  which  the  effluent  concentration  varies  irregularly  and 
tends  to  be  smaller  than  that  obtained  by  the  deep  bed  filter.   This  exception 
was  surprising  in  that  the  shortest  fixed  bed  (25%  of  the  entire  bed  height) 
tended  to  remove  the  largest  amount  of  suspended  solids.   However,  it  was  observed 
that  during  this  run,  the  initial  fixed  bed  was  unstable;  some  of  the  sand  deta- 
ched from  the  bottom  of  the  fixed  bed  into  the  fluidized  portion  of  the  bed,  and 
as  compensation,  some  of  the  fluidized  sand  deposited  and  filled  these  cavities. 
This  detachment-deposition  phenomenon  continued  throughout  the  filtration  run. 
It  is  conceivable  that  this  phenomenon  induces  a  turbulent  retarding  action  on 
the  suspended  solids,  thus  improving  the  removal  efficiency.   No  phenomenon  of 
"cake"  formation  was  observed,  and  the  operation  of  the  semifluidized-bed  filters 
did  not  show  distinct  advantage  over  that  of  the  deep  bed  filter. 

The  suspension  in  Runs  G-1  through  G-3  contained  -  /^250  to  +  //270  (250x270) 
mesh  coal  particles  with  an  average  diameter  of  approximately  0.058  mm.   The 
influent  concentration  was  kept  at  0.0164  wt%  for  each  run.   The  results  of  this 
set  of  runs,  as  shown  in  Figures  VII-5  and  VII-6,  are  quite  similar  to  those 
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obtained  in  Runs  F-1  through  F-3.  Again,  for  the  particles  of  this  size  no 
"cake"  formation  was  observed,  and  the  operation  of  the  semifluldized-bed 
filters  did  not  show  distinct  advantage  over  that  of  the  deep  bed  filter. 

Filtration  of  Relatively  Large  Particles 

In  Runs  H-1  through  H-3,  the  suspension  comprised  -  #140  to  +  #200 

(140x200)  mesh  coal  particles  with  an  average  diameter  of  approximately  0.090 

nm.  The  Influent  concentration  was  maintained  at  0.0149  wt%  for  each  run. 

As  indicated  in  Chapter  VI,  around  99%  of  the  suspended  solids  with  this  size 

was  removed  with  the  5  cm  deep  bed  filter  (see  Figure  VI-10) ,  inferring  that 

the  shallow  penetration  of  suspended  solids  will  be  significant  in  the  operation 

of  the  20  cm  deep  bed  filter.  Filtration  runs  were  performed  with  the  deep  bed 

filter  (Run  H-1) ,  the  semlf luidlzed-bed  filter  with  h  /h  =  25%  (Run  H-2) ,  and 

pa 

the  semifluidlzed-bed  filter  with  h  /h  =  50%  (Run  H-3) .   The  change  in  the 

pa 

pressure  drop  for  each  run  is  plotted  in  Figure  VII-7.   In  operating  the  deep 

bed  filter  (Run  H-1),  a  sharp  Increase  in  the  pressure  drop  was  observed  around 

5  minutes  after  the  onset  of  filtration.   This  is  obviously  Indicative  of  the 

rapid  clogging  around  the  entrance  section  of  the  filter  bed,  that  is,  "cake" 

formation  becoming  significant.  The  results  plotted  in  Figure  VII-7  show  that 

the  capacity  of  the  semifluidlzed-bed  filters,  reflected  in  the  duration  of 

filtration  runs,  was  substantially  higher  than  that  of  the  deep  bed  filter;  the 

duration  of  the  semifluidlzed-bed  filter  with  h  /h  =  25%  (Run  H-2)  was  2  times 

pa 

that  of  the  deep  bed  filter  (Run  H-1) ,  and  a  duration  3  times  that  of  the  deep 

bed  filter  was  attained  with  the  semifluidlzed-bed  filter  with  h  /h  =  50%  (Run 

pa 

H-3),  All  the  filters  achieved  essentially  100%  removal  of  suspended  solids 

except  the  semifluidlzed-bed  filter  with  h  /h  =  25%  which  yielded  around  99% 

pa 

removal  of  suspended  solids.  This  result  is  reasonable  in  that  the  thinner  the 
fixed  section,  the  lower  the  removal  efficiency.  Note,  however,  that  this 
result  is  different  from  that  observed  in  runs  of  series  F;  in  these  runs,  the 
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semifluidlzed-bed  filter  with  the  shortest  fixed  section  tended  to  remove 

the  largest  amount  of  suspended  solids  because  of  the  detachment-deposition 

phenomenon  discussed  previously.   It  appears  that  the  semlf luidized-bed  filter 

with  h  /h  =  50%  had  the  best  performance;  the  pressure  drop  buildup  was  the 

smallest  without  deterioration  in  the  quility  of  the  filtrate. 

The  suspension  in  Run  I-l  through  1-3  consisted  of  -  #100  to  +  #140 

(100x140)  mesh  coal  particles  with  an  average  diameter  of  approximately  0,127  mm. 

The  influent  concentration  was  kept  at  0.0140  wt%  for  each  run.   The  results 

shown  in  Figure  VII-8  resemble  those  in  Figure  VII-7.   For  the  run  with  the  deep 

bed  filter  (Run  I-l) ,  a  sharp  increase  in  the  pressure  drop  was  observed  around 

3  minutes  after  onset  of  the  run.  All  the  filters  achieved  essentially  100% 

removal  of  suspended  solids.   The  duration  of  the  semifluldized-bed  filter  with 

h  /h  =  25%  (Run  1-2)  was  3  times  that  of  the  deep  bed  filter  (Run  I-l),  and  a 
pa 

duration  5  times  that  of  the  deep  bed  filter  was  attained  with  the  semif luidized- 

bed  filter  with  h  /h  =  50%  (Run  1-3) ,   Again,  the  semlf luldlzed-bed  filter  with 
pa 

h  /h  =  50%  exhibited  the  best  performance, 
pa 

Effect  of  the  Size  of  Suspended  Solids 

Note  that  the  size  ratio  of  suspended  solids  to  sand  is  an  important  para- 
meter affecting  the  performance  of  a  semlf luldlzed-bed  filter.  As  discussed  in 
the  previous  sections,  the  larger  the  ratio,  the  more  significant  the  shallow 
penetration  and  "cake"  formation  of  suspended  solids.   It  can  concluded,  there- 
fore, that  when  the  size  ratio  exceeds  a  certain  value,  the  shallow  penetration 
and  "cake"  formation  of  suspended  solids  become  significant,  and  the  use  of  the 
semifluidized-bed  filters  can  substantially  Improve  the  efficiency  of  the  opera- 
tion as  illustrated  in  Figures  VII-7  and  VII-8, 

Structural  Evolution  of  the  Semifluidized-bed  Filters 

The  photographic  study,  presented  in  Figures  VII-9  through  VII-14,  depict 
clearly  the  structural  evolution  of  the  semlfluldlzed-bed  filters  during  opera- 
tion. 


Semlfluidlzed-bed  filter  with  h  /h  =  25%.      Prior  to  a  filtration  run, 
_E2 

the  semifluidized-bed  filter  with  h  /h  =  25%  had  the  configuration  as  seen 

pa  ° 

in  Figure  VII-9.   At  the  onset  of  operation,  a  fixed  section  and  a  fluidized 
section  were  formed  in  series  in  the  semif luidized  bed,  in  which  the  fixed-bed 
height  was  approximately  25%  of  the  entire  bed  height.   Initially,  the  entire 
bed  acted  like  a  deep  bed  filter  with  the  fluidized  section  exerting  relatively 
little  effect  on  the  operation.   As  the  clogging  in  the  fixed  section  became 
appreciable,  the  suspended  solids  could  no  longer  penetrate  into  the  fixed 
section.   With  the  accumulation  of  the  suspended  solids  in  the  fluidized  section, 
the  porosity  of  the  fluidized  section  was  reduced,  rendering  the  linear  velocity 
of  flow  in  this  section  to  increase.   This,  in  turn,  caused  an  increase  in  the 
drag  forces  exerted  on  the  fluidized  particles,  forcing  the  suspended  solids 
and  fluidized  sand  to  form  an  extended  fixed  section  because  the  fluidized 
section  could  not  expand.   Notice  that  a  predominant  particles  circulation 
pattern  is  often  observed  during  the  operation  of  a  fluidized  bed  (see,  e.g., 
Kunil  and  Levensipiel,  1977).   Furthermore,  this  circulation  pattern  often  skews 
to  one  side  of  the  fluidized  bed  because  of  the  difficulty  in  maintaining  the 
completely  uniform  flow  condition  within  the  bed;  the  flow  through  the  distri- 
butor plate  is  seldom  totally  uniformly  distributed  or  the  suspended  solids  are 
not  uniformly  captured  in  the  fixed  section.   This  nonsymmetrical  flow  pattern 
in  the  fluidized  section  resulted  in  a  higher  linear  velocity  of  flow  in  one 
side  of  this  section  [see  Figure  Vll-lO(a)].   The  flow  with  a  higher  linear 
velocity  delivered  an  additional  amount  of  particles  to  the  interface  between 
the  fixed  and  fluidized  sections,  thereby  extending  the  fixed  section  along 
this  side  [see  Figure  Vll-lO(a) ,  10(b)].   This  Increased  the  resistance  to  the 
flow  and  reduced  the  flow  passing  through  it  [see  Figure  Vll-lO(b)],  which, 
in  turn,  altered  the  flow  path,  thereby  causing  the  local  linear  velocity  of 
flow  along  the  extended  fixed  section  to  increase  further  and  promoting  the 
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continuous  extension  of  the  fixed  section  [see  Figure  Vll-lO(b),  10(c)]. 
This  process  continued  during  the  run  until  the  fixed  section  extended  to 
the  bottom  distributor  plate.   It  appears  that  instead  of  forming  a  layer  of 
"cake"  in  the  fixed  section,  most  of  the  suspended  solids  distributes  throughout 
the  extended  fixed  section,  thus  reducing  the  rate  of  the  pressure  drop  increase. 
The  structural  evolution  of  the  semif luidized  bed  during  the  filtration  run 
(Run  1-2)  is  shown  in  Figure  VII-11.   These  phenomena  are  similar  to  those 
observed  by  Hsu  and  Fan  (1981). 

Semifluidized-bed  filter  with  h  /h  =  50%.   Prior  to  a  filtration  run,  the 

pa 

semifluidized-bed  filter  with  h  /h  =  50%  had  the  configuration  as  seen  in  Figure 
pa 

VII-12  or  Figure  VII-13(a).   The  initial  stage  of  the  run  was  similar  to  that 
observed  in  the  operation  of  the  semifluidized-bed  filter  with  h  /h  =  25%. 
However,  instead  of  extending  gradually,  the  fixed  section  extended  rather 
rapidly  along  one  side  down  to  the  bottom  distributor  plate  [see  Figure  VII- 
13(b)].   Because  of  the  high  resistance  to  the  flow  on  this  side  of  the 
completely  fixed  section,  almost  all  the  flow  began  to  pass  through  the  fluidized 
section.   Because  of  this  flow  bypassing,  the  extended  fixed  section  barely 
contributed  to  the  pressure  drop.   As  the  filtration  run  continued,  an  increasing 
amount  of  suspended  solids  and  fluidized  sand  deposited  themselves  on  the  inter- 
face between  the  extended  fixed  section  and  fluidized  section;  subsequently, 
this  action  thickened  the  extended  fixed  section  along  the  horizontal  direction 
[see  Figure  VII-13(c)].   This  lateral  growth  of  the  extended  fixed  section  seemed 
to  be  due  to  the  uneven  flow  motion  in  the  fluidized  section  as  elaborated  in 
the  following.   Referring  to  Figure  VII-13(c),  the  left  side  of  the  fixed  section 
was  highly  resistant  to  the  flow  because  of  a  large  amount  of  deposition  of 
suspended  solids.   Furthermore,  the  right  side  of  the  fixed  section  was  also 
highly  resistant  to  the  flow  because  of  the  formation  of  the  completely  fixed 
section.   Therefore,  the  flow  in  the  fluidized  section  tended  to  pass  through 
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the  corner  formed  by  the  initial  fixed  section  and  the  extended  fixed  section 
as  illustrated  in  Figure  VII-13(c);  this  was  the  direction  the  least  resistant 
to  the  flow.   The  lateral  component  of  this  flow  would  constantly  deliver  both 
the  suspended  solids  and  fluidized  sand  to  the  interface  between  the  extended 
fixed  section  anf  the  fluidized  section,  thus  resulting  in  the  continuous  growth 
of  the  extended  fixed  section.   It  appears  that  this  lateral  growth  of  the 
extended  fixed  section  captures  continuously  both  the  suspended  solids  and 
fluidized  sand;  simultaneously,  the  flow  tends  to  bypass  this  growing  fixed 
section,  thus  preventing  a  significant  increase  in  the  pressure  drop.   A 
structural  evolution  of  the  semif luidized  bed  during  a  filtration  run  (Run  1-3) 
is  illustrated  in  Figure  VII-14, 
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Figure  VII-2.      Photograph  of  the  semif luidized-bed    filter 
during  the  course  of  backwashlng. 
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Figure  VII-3.      Increase    in   the   pressure   drop  as   a    function  of    filtration   ti™    : 
suspension   concentration,    0.0227  wt/.,  ,n   R  V  ■ 

Runs;    F-1,    entire   bed   height,    20  cm;    pH ,    8.0;    T.    20   8  C 

F-2      entire   bed   height,    23.6   cm;    pH,    8.1;    T,    21. UO, 
F-3'    entire   bed   height,    22.4   cm;    pH ,    8.0;    T,    20.7   C. 


VII-16 


30 


25 


20 


c 
o 

2|5 

■♦- 
c 
<u 
o 
c 
o 
o ,  ^ 


UJ 


O         h     /h  =  1007.    (Deep  Bed)  :    Run  F-1 
pa 

O        h     /h  =  50y.:    Run  F-3 
pa 

A         h     /h  =  257.:    Run  F-2 
pa 


oo 


o 


o  o  o  o 


o 


o 


o 


§ 


o 


a  <^ 


6 


I    I    1 


_L_L 


12 


14 


18        20 


Filtration  Time   (min) 


flmre  VII-4        Effluent   concentration  of   suspended    particles   as   a    function 
^  of    filtration    time    ;    suspension  concentration,    0.0227  wt/.; 

Runs   F-1   through  F-3. 
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Figure  VII-5.      Increase    in    the    pressure   drop  as    a    function  of    filtration    time 
suspension  concentration,    0.0164  wt/.;  ^ 

Runs;   G-1,   entire   bed  height,    20  cm;    pH ,   8.5;   T,    22.1   C ;_ 

G-2,   entire  bed   height,    23.6  cm;    pH ,   8.5;   T,    21.9  _C; 

G-3,    entire    bed   height,    22.4   cm;    pH ,    8.5;    T,    21.8  C. 
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Flaure  VII-6       Effluent   concentration  of  suspended   particles  as   a   function 
'      of   filtration   time    :    suspension  concentration,   0.0164  wt7.; 
Runs  G-1   through  G-3. 
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neure  VII-7.      Increase    in    the    pressure   drop  as   a    function   of    filtration    time 
suspension  concentration,    0.0149  wt%; 
Runs;    H-1,    entire  bed   height,    20  cm;    pll ,    9.4;   T,    16.6  C -^ 

H-2,   entire   bed   height,    23.6   cm;    pH ,    9.3;   T,    16.5_C; 

H-3,    entire   bed   height,    22. A    cm;    pH ,    9.3;    T,    16.4   C. 


VII- 20 


400 
350 

300 


CD 
-250 


O 

E  200 

o 

Q. 

o 

^— 
Q 


150 


=3 
<V 


1-2 


"       oo 


6 


O 


Ao 


A     o^ 


5: 


o  h      /h  =  100%    (Deep  Bed) :    Run  I-l 

'^  pa 

A  h     /h  =  50%:    Run  1-3 

^/  pa 

^  h      /h  =  25%:    Run  1-2 


100 


80 


^ 


J_ 


_L 


J 


0     2 


8      10      12      14       16     18      20    22     24     26    28     30 

Filtration  Time    (min) 


Figure  VII-8. 


Increase    in    Che    pressure  drop  as   a    function   of    filtration    time 

suspension  concentration,    0.0140  wt7. ; 

Runs;    I-l,    entire   bed   height,    20   cm;    pH ,    9.1;    T,    17.0  C ;_ 

1-2,    entire   bed   height,    23.6   cm;    pH ,    9.0;    T,    17.1   C; 

1-3      entire   bed    height,    22.4   cm;    pH ,    9.1;    T,    17.2   C. 
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Figure  VII-9.   Photograph  of  the  semlf luldized-bed  filter  prior 

to  the  start  of  a  filtration  run:  h  /h  "  257.. 
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Figure  VII-12.   Photograph  of  the  semlf luidized-bed  filter  prior 

to  the  start  of  a  filtration  run:  h  /h  =  50%. 
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CHAPTER  VIII   CONCLUSIONS  AND  RECOMI-IENDATIONS 

In  this  chapter,  significant  conclusions  of  this  work  are  sunmiarized 
and  possible  future  works  are  recommended, 

CONCLUSIONS 

The  deep  bed  filtration  is  an  extremely  complex  phenomenon  involving 
numerous  linear  and  nonlinear  mechanisms.   For  such  a  complex  system,  it 
is  impossible  to  construct  an  all-inclusive  model  which  is  capable  of  taking 
into  account  all  known  nonlinear  and  complex  mechanisms.   In  fact,  any  model 
which  attempts  to  describe  every  mechanism  becomes  excessively  complex  and 
useless.   The  basic  tenet  of  the  modeling  is  that  the  resultant  model  should 
be  as  simple  as  possible  and  yet  sufficiently  precise  and  mechanistic  to 
represent  or  correlate  the  experimental  data  within  their  accuracy.   In 
other  words,  to  derive  a  model  for  such  a  highly  complex  or  nonlinear  system, 
we  can  not  but  resort  to  the  gross  or  linear  approximation  to  some  degree. 
Needless  to  say,  a  model  of  any  system  is  almost  always  operational  in  that 
it  should  be  concerned  with  only  the  limited  number  of  the  significant 
parameters  or  state  variables  of  the  system,  which  are  of  interest  to  the 
modeler. 

In  this  work,  the  pressure  drop  characteristics  through  the  deep  bed 
are  of  primary  concern.   To  simulate  the  pressure  drop  dynamics  under  a 
constant  flow  condition,  stochastic  models,  namely,  the  linear,  birth-death 
process  and  second-order,  pure  birth  process,  in  conjunction  with  the  Carman- 
Kozeny  equation,  have  been  proposed.   These  models  take  into  account  the 
blockage  of  pores  by  suspended  solids  and/or  scouring  of  deposited  solids ; 
the  filter  bed  is  assumed  to  be  spatially  lumped  or  to  be  composed  of  several 
compartments.   In  spite  of  this  simplicity,  the  present  models  appear  to 
represent  the  majority  of  the  experimental  data. 


The  resultant  model  equations,  Eqs.  IV-28  and  V-43,  are  simple  and 
easy  to  manipulate.   Three  parameters,  a  and  S  in  Eq.  IV-28  and  k  In  Eq. 
V-43,  are  easily  Identifiable.   By  resorting  to  these  parameters,  the 
influence  of  various  operating  variables  on  the  pressure  drop  characteristics 
can  be  correlated  with  reasonable  ease.   These  variables  include  the  size 
distributions  of  collector  grains  and  suspended  solids,  properties  of  the 
liquid,  filtration  rate,  bed  characteristics,  and  suspension  concentration. 
For  design  and  scale-up  purposes,  the  establishment  of  relationships  between 
each  of  the  parameters  and  the  operating  variables  is  the  key  that  renders 
the  present  models  useful.   Based  on  this  consideration,  much  effort  has 
been  spent  to  experimentally  Investigate  the  relationships  with  the  size 
distribution  of  suspended  solids  and  suspension  concentration  as  presented 
in  Chapter  VI. 

The  results  of  the  semif luidlzed  bed  filtration  show  that  the  capacity 
of  semifluldized  bed  filters,  expressed  in  terms  of  filtration  run  duration, 
ranged  from  2  to  5  times  that  of  the  deep  bed  filter.   Furthermore,  no 
appreciable  differences  were  observed  in  the  quality  of  effluents  obtained 
through  the  two  types  of  filters.   The  results  also  indicate  that  the 
duration  of  semifluldized  bed  filtration  runs  is  strongly  dependent  on  the 
dynamic  mechanism  occurring  in  the  fluidized  portion  of  the  filter  as  illus- 
trated in  Figures  V-7  through  V-12. 

RECOMMENDATIONS 

The  use  of  the  semifluldized  bed  filter  holds  great  promise  for  Increasing 
"^^   filtration  efficiency.   To  aid  Che  design  and  scale-up  of  such  filters 
through  experimental  observations  and  modeling,  the  following  is  recommended. 
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(1)  To  determine  the  performance  of  the  semlfluidlzed-bed  filter,  a  companion 
study  of  the  deep  bed  filter  has  been  carried  out  to  establish  a  basis 
for  comparison  and  evaluation.   Stochastic  models  have  been  developed  to 
simulate  the  pressure  drop  change  through  the  deep  bed  filter.   However, 
the  usefulness  of  these  models  as  such  in  the  design  and  scale-up  of  the 
deep  bed  filter  is  rather  limited.   For  these  purposes,  the  establishment 
of  relationships  between  each  of  the  parameters,  a,    g,  and  k,  and  the 
operating  variables,  including  the  size  distributions  of  grains  and  sus- 
pended solids,  size  ratio  of  suspended  solids  to  grains,  filtration  rate, 
bed  characteristics,  and  suspension  concentration,  is  the  key  that  renders 
the  present  models  useful.   Thus,  it  is  necessary  to  carry  out  an  extensive 
experimental  study  of  the  deep  bed  filtration  under  a  variety  of  operating 
conditions  to  determine  empirically  or  semiempirically  such  relationships. 

(2)  The  present  models  do  not  take  into  account  the  variation  in  the  specific 
deposit  of  suspended  solids.   The  suitability  of  neglecting  the  effect  of 
this  variation  need  be  further  justified.   Furthermore,  the  present  models 
are  unable  to  simulate  the  phenomenon  of  cake  formation  dominent  in  some 
situations;  additional  modeling  work  should  be  carried  out  to  complement 
this  deficiency, 

(3)  The  results  of  semlfluidlzed-bed  filtration  show  that  the  capacity  of  the 
semifluidlzed-bed  filter  can  be  Increased  up  to  several  folds  that  of  an 
equivalent  deep  bed  filter.   This  need  not  be  the  limit;  for  example, 
this  improvement  will  be  even  greater  by  augmenting  the  depth  of  deep  bed 
employed  in  this  work,  which  is  only  20  cm.   To  determine  the  feasible 
and/or  optimal  regions  of  operation,  an  experimental  study  on  the 
semlfluidlzed-bed  filter  should  be  carried  out  by  systematically  varying 
the  operating  variables,  including  the  entire  bed  height,  degree  of  flui- 
dization,  filtration  rate,  size  distributions  of  grains  and  suspended  solids 


and  size  ratio  of  suspended  solids  to  f^ralns.  *  -^  "" 

(4)  A  rather  high  flow  rate,  greater  than  13.1  m/h,  was  necessary  In  malntal-  . 
nlng  the  fluldlzatlon  of  grains  In  the  semlf luldlzed-bed  filters  In  this 
work.   On  the  other  hand,  for  the  purification  of  potable  water  and  the 
treatment  of  waste  water,  filtration  runs  are  usually  performed  under  the 
flow  rate  In  the  neighborhood  of  5  m/h  and  with  the  suspended  solids  of 
the  average  diameter  between  1  and  50  ym.   To  operate  semlf luldlzed-bed 
filters  under  such  conditions,  particles  smaller  and  lighter  than  the 
20x30  mesh  sand  used  in  this  work  ought  to  be  employed  as  the  filtering 
medium. 

(5)  The  suspended  solids  In  the  potable  water  or  waste  water  often  vary  widely 
in  size  either  by  the  flocculatlon  from  the  secondary  effluent  or  by  the 
addition  of  filter  aid.   Recall  that  the  pressure  drop  Increase  was  quite 
large  for  the  deep  bed  filtration  with  suspended  solids  varying  widely 

in  size  because  the  larger  solids  tended  to  clog  the  bed,  preventing  the 
penetration  of  smaller  solids  into  the  bed.   However,  it  may  be  possible 
that  in  filtering  such  suspended  solids  with  the  semlfluldlzed  bed,  the 
larger  solids  will  tend  to  clog  the  bed  and  cause  the  extension  of  the 
fixed  section;  this  extended  fixed  section  will  aid  substantially  the 
capture  of  the  smaller  solids,  reducing  the  rate  of  pressure  drop  increase. 
Therefore,  it  is  possible  that  the  improvement  by  the  semlf luldlzed-bed 
filter  will  be  even  greater  than  that  observed  for  the  filtratlons  with 
suspended  solids  with  a  narrow  size  range.   This  possibility  should  be 
experimentally  demonstrated. 

(6)  A  mathematical  model  of  the  semlf luldlzed-bed  filtration  should  be  developed 
for  the  design  and  scale-up  of  the  filter.   At  the  initial  stage  of  filtra- 
tion, the  semlfluldlzed  bed  will  behave  like  a  deep  bed;  the  present  models 
for  the  deep  bed  filtration  are  capable  of  simulating  its  performance  at 
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this  stage.  As  the  clogging  In  the  fixed  section  becomes  appreciable,  the 
fixed  section  will  start  extending.  In  modeling  this  stage  of  operation, 
tvo  parameters  are  important,  one  being  the  rate  of  the  dovmward  extension 
of  fixed  section  and  the  other  being  the  rate  at  which  the  suspended  solids 
are  deposited  through  this  extended  fixed  section.  The  estimation  of  these 
rates  may  result  in  a  model  for  depicting  the  pressure  drop  increase  during 
this  stage.   Eventually,  the  fixed  section  will  extend  down  to  the  bottom 
of  the  bed,  signaling  the  onset  of  the  last  stage  of  operation.  At  this 
stage,  the  lateral  component  of  uneven  flow  in  the  fluldlzed  section,  as 
discussed  in  Chapter  VII,  will  force  suspended  solids  and  fluldlzed  sand  to 
be  captured  by  the  extended  fixed  section  and  thicken  it  along  the  horizontal 
direction.  Simultaneously,  the  flow  will  tend  to  bypass  this  extended  fixed 
section  because  of  its  high  resistance  to  the  flow,  Indicating  that  this 
section  will  contribute  little  to  the  pressure  drop  increase.   Most  of  the 
pressure  drop  increase  can  be  attributed  to  the  suspended  solids,  not  captured 
by  the  extended  fixed  section,  but  deposited  mainly  in  the  initial  fixed 
section  through  which  the  flow  pass. 
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The  present  study  aimed  to  carry  out  a  comprehensive  investigation 
of  the  performance  of  the  semif luidized-bed  filter.   Because  the  filtration 
In  a  semifluidized-bed  Involves  the  deep  bed  filtration,  a  reasonable  model 
of  semlfluldized-bed  filtration  requires  the  incorporation  of  a  deep  bed 
filtration  model.   Therefore,  a  companion  study  on  the  deep  bed  filter  has 
also  been  carried  out. 

Stochastic  models,  namely,  the  linear,  birth-death  process  and  the 
second-order,  pure  birth  process,  have  been  employed  in  conjunction  with  the 
Carman-Kozeny  equation  to  simulate  the  performance  of  the  deep  bed  filter  in 
terms  of  the  pressure  drop  dynamics  under  a  constant  flow  condition.   These 
models  take  into  account  the  blockage  of  pores  by  suspended  solids  and/or 
scouring  of  deposited  solids;  the  filter  bed  is  assumed  to  be  spatially 
lumped  or  to  be  composed  of  several  compartments.   The  present  models  appear 
to  represent  the  majority  of  the  experimental  data. 

A  series  of  experimental  runs  of  deep  bed  filtration  was  conducted 
with  suspensions  of  solids  of  various  sizes,  size  distributions  and  concen- 
trations.  For  each  run,  the  pressure  drop  and  effluent  concentration  were 
measured  as  functions  of  filtration  time.   The  results  are  treated  and  ana- 
lyzed stochastically. 

The  operation  of  deep  bed  filters  is  often  inefficient;  it  is  typically 
characterized  by  the  shallow  penetration  of  suspended  solids  and  rapid  clogging 
in  the  bed.   Fortunately,  these  disadvantages  can  be  overcome  by  the  use  of 
semifluidized-bed  filters.   The  results  of  an  experimental  study  show  that  the 
filtration  capacity  of  the  semifluidized-bed  filter  is  considerably  greater 
than  that  of  an  equivalent  deep  bed  filter.   In  addition  to  pressure  drop  and 
effluent  quality  data,  a  series  of  photographs  is  Included  to  describe  and 
explain  the  mechanisms  of  the  semifluidized-bed  filter,  giving  rise  to  its 
superior  performance  over  the  deep  bed  filter. 
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measured  as  functions  of  filtration  time.   The  results  are  treated  and  ana- 
lyzed stochastically. 

The  operation  of  deep  bed  filters  is  often  inefficient;  it  is  typically 
characterized  by  the  shallow  penetration  of  suspended  solids  and  rapid  clogging 
in  the  bed.   Fortunately,  these  disadvantages  can  be  overcome  by  the  use  of 
semifluidized-bed  filters.   The  results  of  an  experimental  study  show  that  the 
filtration  capacity  of  the  semifluidized-bed  filter  is  considerably  greater 
than  that  of  an  equivalent  deep  bed  filter.   In  addition  to  pressure  drop  and 
effluent  quality  data,  a  series  of  photographs  is  Included  to  describe  and 
explain  the  mechanisms  of  the  semifluidized-bed  filter,  giving  rise  to  its 
superior  performance  over  the  deep  bed  filter. 


